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Abstract. Shear viscosity is a measure of the amount of dissipation in a simple fluid. 
In kinetic theory shear viscosity is related to the rate of momentum transport by quasi- 
particles, and the uncertainty relation suggests that the ratio of shear viscosity rj to 
entropy density s in units of h/ks is bounded by a constant. Here, h is Planck's 
constant and ks is Boltzmann's constant. A specific bound has been proposed on 
the basis of string theory where, for a large class of theories, one can show that 
■q/s > U/^Trks)- We will refer to a fluid that saturates the string theory bound as a 
perfect fluid. In this review we summarize theoretical and experimental information 
on the properties of the three main classes of quantum fluids that are known to have 
values ofrj/s that are smaller than h/ks- These fluids are strongly coupled Bose fluids, 
in particular liquid helium, strongly correlated ultracold Fermi gases, and the quark 
gluon plasma. We discuss the main theoretical approaches to transport properties of 
these fluids: kinetic theory, numerical simulations based on linear response theory, and 
holographic dualities. We also summarize the experimental situation, in particular 
with regard to the observation of hydrodynamic behavior in ultracold Fermi gases and 
the quark gluon plasma. 
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1. Introduction 

A fluid is a material that can be described by the laws of fluid dynamics. These laws 
imply that the response of a fluid to slowly varying external perturbations is completely 
governed by conservation laws. In the case of simple fluids, such as water, the conserved 
quantities are mass, energy, and momentum. 

The study of fluids is one of the oldest problems in physics p. Understanding why 
certain materials make good fluids, and others do not, has nevertheless remained a very 
difficult question. The quality of a fluid can be characterized by its shear viscosity rj. 
Shear viscosity is defined in terms of the friction force F per unit area A created by a 
shear flow with transverse flow gradient V y v x , 



Viscosity causes dissipation which converts part of the kinetic energy of the flow to heat. 
A good fluid is therefore characterized by a small shear viscosity. Indeed, the inverse of 
shear viscosity, ip = 1/r], is sometimes called fluidity. 

The molecular theory of transport phenomena in dilute gases goes back to Maxwell. 
Maxwell realized that shear viscosity is related to momentum transport by individual 
molecules. A simple estimate of the shear viscosity of a dilute gas is 



where n is the density, p is the average momentum of the molecules, and l m f p is the 
mean free path. The mean free path can be written as l m f p = l/(na) where a is a 
suitable transport cross section. This implies that the shear viscosity of a dilute gas 
grows with temperature (as p ~ T 1 / 2 ) but is approximately independent of density. This 
counterintuitive result is confirmed by experiment, going back to experiments carried 
out by Maxwell himself [2]. Equ. ([2]) also shows that the viscosity of an ideal gas is 
infinite, not zero. In order to achieve thermal equilibrium we have to view the ideal gas 
as the limit of an interacting system in which the scattering cross section a is taken to 
zero. In this limit the mean free path, and with it the viscosity, goes to infinity. 

At low temperature gases condense into the liquid (or solid) state. In a liquid 
transport is no longer governed by the motion of individual molecules. A simple picture, 
due to Frenkel, Eyring, and others, is that momentum transport is due to processes that 
involve the motion of vacancies [3] . These processes can be viewed as thermally activated 
transitions in which a molecule or a cluster moves from one local energy minimum to 
another. The viscosity scales as [U [5] 



— = T] VyV X . 



(2) 



7] ~ hne E/{ksT) , 



(3) 
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Table 1. Viscosity ?y, viscosity over density, and viscosity over entropy density ratio 
for several fluids. Data for water and helium taken from [BJ [7] and [8], data for Li and 
the quark gluon plasma (QGP) will be explained in Sect.[5l For water and helium we 
show data at atmospheric pressure and temperatures just below the boiling point and 
the A transition, respectively. These data points roughly correspond to the minimum 
of rj/n at atmospheric pressure. We also show and data near the tri-critical point which 
roughly corresponds to the global minimum of rj/s. Note that the quark gluon plasma 
does not have a well defined density. 

where E is an activation energy and h is Planck's constant. We note that the viscosity of 
a liquid has a very strong dependence on temperature. We also observe that the overall 
scale involves Planck's constant. The appearance of h is related to Eyring's assumption 
that the collision time of the molecules is h/(kBT), the shortest time scale in a liquid. 
We will come back to this assumption below. Equ. (EJ) shows that the viscosity of a 
liquid grows as the temperature is lowered. Together with equ. (jlj) this result implies 
that the viscosity of a typical fluid has a minimum as a function of temperature, and 
that the minimum is likely to occur in the vicinity of the liquid-gas phase transition. 

Experimental results show that the minimum value of the viscosity of good fluids, 
like water, liquid helium and liquid nitrogen, differs by many orders of magnitude, see the 
data in Table [TJ The SI unit for viscosity is Pascal second (Pa • s), the CGS unit is Poise 
(P). Note that 1 Pa • s = 10 Poise. Clearly, it is desirable to normalize the viscosity to a 
suitable thermodynamic quantity in order to make more useful comparisons. Equations 
(T2]) and indicate that a suitable ratio is provided by rj/n. We note that the ratio of 
viscosity over mass density p = mn is known as the kinematic viscosity v — rj/p. The 
behavior of solutions of the Navier-Stokes equation is governed by the Reynolds number 

Re = (£j vL , (4) 

which is the ratio of a property of the flow, its characteristic velocity v multiplied by its 
characteristic length scale L, over a property of the fluid, its kinematic viscosity. Good 
fluids attain larger Reynolds numbers, and are more likely to exhibit turbulent flow. 
Data for the ratio r)/n are tabulated in Table [TJ We observe that the ratios r)/n for 
good fluids are indeed similar in magnitude. 

A disadvantage of considering the ratio r)/n is that it is not possible to include 
relativistic fluids in the comparison. In the case of a relativistic fluid the number of 
particles is not conserved. As a consequence the quantity n is not be well defined in 
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an interacting system. In a quark gluon plasma, for example, only the net number of 
quarks (the number of quarks minus the number of anti-quarks) is well defined, but 
the number of quarks or the number of gluons is not. In Sect. 13.11 we will show that 
the Reynolds number of a relativistic fluid is defined in terms of the ratio n/(sT), 
where s is the entropy density and T is the temperature. This indicates that we should 
consider the ratio rj/s instead oirj/n. We note that this ratio is well defined in both the 
relativistic and non-relativistic limit, and that s ~ nks for many fluids. For example, 
in a non-interacting relativistic Bose gas s/n ~ 3.6 ks, and for non-interacting fermions 
s/n ~ 4.2 ks- In a weakly interacting non-relativistic gas s/n ~ ks up to logarithms of 
gn/(mT) 3 / 2 , where g is the degeneracy factor. Data for n/s in units of h/ks are also 
given in Table HJ 

We observe that good fluids are characterized by n/s ~ %/ks- This value is 
consistent with simple theoretical estimates. Consider the kinetic theory estimate in 
equ. ([2]). In the strong coupling limit the mean free path becomes very small, but the 
uncertainty relation suggests that pl m f p ^ Ti [9]. For a rough estimate we may use the 
relation s ~ 3.6ksn and conclude that r]/s > h/(10.8 k^). A bound on n/s is also 
indicated by the high temperature limit of Eyring's formula equ. (131). We note, however, 
that kinetic theory is not reliable in the regime n/s ~ h/ks, and other methods are 
needed to determine the minimum value of n/s. 

A precise value of the viscosity bound was proposed based on results from string 
theory. Kovtun et al. conjectured that [10] 

for all fluids. We will call a fluid that saturates the bound ([5]) a "perfect fluid" . A perfect 
fluid dissipates the smallest possible amount of energy, and satisfies the laws of fluid 
dynamics in the largest possible domain. In a typical fluid, hydrodynamics is an effective 
description of long wavelength fluctuations, but in a perfect fluid hydrodynamics is 
reliable at distances as short as the inter-particle spacing. 

The viscosity bound conjecture raises a number of interesting questions: 

• Is the conjecture in equ. ([5]) correct? Does this bound, or some other bound on r]/s, 
follow from the general principles of quantum mechanics and statistical mechanics? 

• Is there a "perfect fluid" in nature, i.e. can we observe a fluid that attains the value 
n/s = ^/(47tA;b)? If yes, what are the characteristics of such a fluid? Is it possible 
to describe the fluid in terms of quasi-particles? 

• How is n/s correlated with other transport properties, like bulk viscosity, diffusion 
constants, conductivities? Are there bounds on other transport properties? 

We will not be able to provide definitive answers to these questions in this review. 
There are, however, a number of recent results, from both theory and experiment, that 
shed light on these issues: 
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• The experimental realization of new classes of quantum fluids. Prior to 1995 the 
only bulk quantum fluids that could be studied in the laboratory were the two 
isotopes of liquid helium, 4 He and 3 He. In 1995 several groups achieved quantum 
degeneracy is dilute atomic Bose gases. In 1999 experimentalists also succeeded in 
producing degenerate atomic Fermi gases pj], [12] • Using Feshbach resonances it is 
possible to experimentally control the interaction between the atoms, and to study 
equilibrium and transport properties as a function of the interaction strength. 

• The experimental discovery of almost ideal hydrodynamic flow in a completely 
different physical system, the quark gluon plasma created in heavy ion collisions 
at the relativistic heavy ion collider (RHIC) at Brookhaven National Laboratory 
[TBI [TU [15]. The quark gluon plasma also exhibits a large energy loss for high 
energy colored particles, and a very small heavy quark diffusion constant. 

• Progress in non-equilibrium field theory culminated in the calculation of transport 
coefficients of weakly coupled gauge theory plasmas [HI [T7J [IB] . These results 
complete the program of using kinetic theory to calculate the transport properties 
of the three main classes of quantum liquids: Bose gases, Fermi gases, and gauge 
theory plasmas. 

• The theoretical discovery of a completely new method for computing the transport 
properties of very strongly coupled fluids [19]. This method is based on the 
holographic duality between certain strongly coupled field theories in d = 4 space- 
time dimensions and weakly coupled string theory in d = 10 [20]. For gauge theories 
that have a weakly coupled string dual the shear viscosity to entropy density ratio 
at infinite coupling is rj/s = ^/(47tA;b). It was also shown that the first correction 
to this result at finite but large coupling increases rj/s [2T] , and it was conjectured 
that rj/s = h/(47rkB) is a universal lower bound [10J. 

It is the goal of this review to summarize these recent developments. For this 
purpose we shall concentrate on three representative fluids: 4 He, a strongly coupled 
Bose fluid; atomic Fermi gases near a Feshbach resonance, which are the most strongly 
coupled Fermi liquids; and the quark gluon plasma near the critical temperature for 
condensation into hadron gas, which is a very strongly coupled plasma. The review 
is structured as follows: In Sect. [2] we discuss the thermodynamics of these quantum 
fluids. In Sect. [3] we review theoretical approaches to transport properties. We briefly 
summarize the hydrodynamic description of relativistic and non-relativistic fluids, as 
well as superfluids in Sects. 13.1113.31 A general connection between transport coefficients 
and the underlying field theory is provided by Kubo relations, which we introduce in 
Sect. 13.41 In Sect. 13.51 we concentrate on fluids that can be described in terms of weakly 
coupled quasi-particles. In this case transport properties can be computed using kinetic 
theory. In Sect.HJwe summarize results for transport coefficients that have been obtained 
using holographic dualities. Finally, in Sect. [5] we discuss experimental results for the 
viscosity and other transport properties of strongly coupled quantum fluids. 

Needless to say, a review of this size cannot adequately summarize all the work 
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that has been done on the transport properties of quantum fluids. A standard reference 
on the properties of liquid helium is [8] , recent reviews on strongly coupled Fermi gases 
are [22l [23] , and the physics of the strongly coupled quark gluon plasma is discussed 
in [23]. The kinetic theory of dilute Bose and Fermi gases is covered in textbooks, 
see [251 26], and the kinetic theory of gauge fields was reviewed in [27]. Reviews of 
the AdS/CFT correspondence with an emphasis on transport theory are [281 129], and 
reviews of relativistic hydrodynamics can be found in [30], [311 [32] . 



2. Strongly coupled quantum fluids 

In this section we will discuss equilibrium properties of strongly interacting quantum 
fluids. We will specify the effective action for bosonic, fermionic, and gauge theory fluids, 
identify the relevant physical scales, and discuss the nature of low energy excitations. 



2.1. Bose fluids: Dilute Bose gases 

A gas of atoms satisfying Bose statistics can be described in terms of a scalar field 
t/>(x, t) governed by the Hamiltonian 

H= Jd 3 x ^*(x, t) t) (6) 

+ \S d3xi I d3x2 ^*( x i>*M x i>*M x i -x 2 )^*(x 2 ,tMx 2 ,t), 

where m is the mass of the boson, and V^(x) is a potential. The Hamiltonian is 
invariant with respect to translations and rotations, as well as under the U(l) symmetry 
ip — > exp(ia)ip. Symmetries correspond to conservation laws. Translations and rotations 
correspond to the conservation of linear and angular momentum, and the U(l) symmetry 
is associated with the conservation of the number of atoms. The Schrodinger equation 
is also invariant under Galilean transformations x — > x — vt which act on the field as 
^(x, t) — > exp(imv • x — |mv 2 t)?/>(x — vt, t). This symmetry will play a role when we 
consider the motion of fluids. 

If the typical momenta are small compared to l/r , where r is the range of 
the potential, we can approximate the interaction by a contact term. For very small 
momenta the leading contribution is an s-wave interaction 

V(x 1 -x 2 ) = C7 5(x 1 -x 2 ), (7) 

where Co can be related to the scattering length a, Co = (AirTi 2 a)/m. In order to 
make connections with the physics of relativistic fluids it is also useful to introduce the 
lagrangian £ = ihtpdoip — Ti, where TC is the Hamiltonian density. The lagrangian is 

c = tf(*m + ^)i>-%(M) 2 . (8) 

In the following we shall consider many-body systems described by this lagrangian. We 
first study the relevant scales in a weakly interacting Bose gas governed by the s-wave 
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scattering length. At high temperature the Bose gas is a classical Boltzmann gas. The 
average energy of the atoms is §&bT, and the average momenta are of order (m/c^T) 1 / 2 . 
The importance of quantum statistics is governed by the parameter uvq, where n is the 
density, vq = A 3 is the quantum volume and 

A - (9) 

is the thermal wave length. Quantum statistics becomes important for uvq ~ 1, and 
Bose condensation in an ideal gas occurs at uvq = 2.61, corresponding to a critical 
temperature 

^^(cOw) ■ (10) 
The effects of a non-zero scattering length can be taken into account order by order in 
an expansion in an 1 / 3 . At high temperature this is the standard virial expansion 

P = nk B T{l + b 2 n + 0(n 2 )} , (11) 

where P is the pressure and b 2 is the second virial coefficient. In the limit of small a 
the second virial coefficient of a single component Bose gas is given by [33] 

b 2 = -^\ 3 + 2a\ 2 , (12) 

where the first term is due to quantum statistics, and the second term is related to 
the interaction. The second virial coefficient is finite in the limit of a large scattering 
length. As a — > oo the interaction part approaches —a/2 A 3 , and the role of interactions 
is governed by the same parameter that controls the effects of quantum statistics. 

The interaction also shifts the critical temperature for Bose condensation. The 
calculation of this shift is a non-perturbative problem, even if the scattering length is 
small. This is related to the fact that fluctuations become large in the vicinity of a 
second order phase transition, and perturbation theory breaks down. One can show 
that 5T C ~ an 1 / 3 ^, and there is a physical argument that a repulsive scattering length 
(a > 0) increases the value of T c [3j]. A numerical calculation using the Landau- 
Ginzburg effective lagrangian gives [35] 

AT C = (1.32 ± 0.02)an 1/3 T c . (13) 

Dilute Bose gases in which the scattering length is attractive are not stable, but it 
is possible to create metastable systems confined by external fields. Weak repulsive 
interactions increase the transition temperature but suppress the condensate fraction. 
If an 1//3 is large then the Bose fluid will typically solidify, but the phase structure 
and critical density depend on details of the interaction. A hard sphere gas freezes 
at an 1 ' 3 ~ 0.24 [36]. 

An issue which is very important for transport properties is the nature of the quasi- 
particle excitations. At high temperature the cross section for binary scattering between 
the atoms decreases with the thermal wavelength, a ~ A 2 ~ l/(mT), and the mean free 
path l m f p ~ l/{na) is large. As a consequence the the atoms are good quasi-particles. 
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At very low temperature the system is superfluid and there is a Goldstone boson, the 
phonon, related to the breaking of the U(l) phase symmetry. Phonons are derivatively 
coupled and the interaction at low energy is weak. This implies that the mean free path 
at low temperature, T <C T c , is also large. The phonon dispersion relation in a weakly 
non-ideal (na 3 < 1) Bose gas was first computed by Bogoliubov. The result is 




For small momenta the dispersion relation is linear, e p ~ c s p, and the phonon velocity 
is given by c s = \/Amanjm. 



2.2. Bose fluids: Liquid helium 

A simple s-wave interaction is sufficient for understanding the properties of trapped 
atomic Bose gases, but more accurate potentials are required for even a qualitative 
description of liquid 4 He. Accurate 4 He potentials can be written as the sum of a short 
range term and a long range van der Waals potential, 

V(r) = V 3r (r)-^. (15) 

The coefficient Cq defines the van der Waals length scale Wdw = (mC^/a 2 ) 1 / 4 . Accurate 
parametrizations of V sr can be found [37J EH] • These potentials have a van der Waals 
length Ivdw — 10.2 do, an effective range r ~ 14 ao, and a very large scattering length 
a ~ 189 ao, where ao = 0.529 A is the Bohr radius. The large s-wave scattering length 
is related to the existence of a very weakly bound 4 He dimer. The binding energy of 
the dimer is B — — 1.1 • 10~ 7 eV. There are many interesting universal effects governed 
by the large scattering length [39]. The density of liquid 4 He is too large for these 
phenomena to be important, but universal effects have been observed in dilute atomic 
gases in which the scattering length is large. 

In the case of 4 He the interaction between the atoms is not weak, and it cannot be 
characterized in terms of the scattering length only. In the high temperature limit 4 He 
is a classical gas, and corrections to the ideal gas behavior are described by the virial 
expansion. The virial expansion provides a very accurate description of the equation of 
state at normal pressure for temperatures above 10 K. At temperatures below 10 K one 
has to rely on quantum Monte Carlo methods or variational many-body wave functions 
[30] . At atmospheric pressure 4 He liquefies at 4.22 K, and it becomes superfluid at 
T c = 2.17 K. This temperature can be compared to the critical temperature for Bose 
condensation of an ideal gas with the density of liquid helium, n = 1/(3.6 A) 3 , which 
is T c ° = 3.1 K. The dependence of T c on the density in the case of a hard sphere gas 
was studied by Griiter et al. [UJ. Helium is well described by an effective hard sphere 
parameter a = 2.20 A. Griiter et al. show that for na 3 < 0.1 the critical temperature is 
larger than that of a non- interacting gas, in agreement with equ. ( TTBl) . The increase in 
T c is small, reaching about 6%. For larger values of na 3 the critical temperature drops 
rapidly, until freezing occurs at na 3 ~ 0.25. 
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The presence of strong interactions also manifests itself in a small condensate 
fraction. Glyde et al. measured the number of condensed atoms N (T) using neutron 
scattering on liquid 4 He at saturated vapor pressure [42J. They find N (T)/N = 
f(l - (T/T c )t with / ~ (7.25 ± 0.75) • 10~ 2 and 7 = 5.5 ± 1. The superfluid transition 
is in the universality class of the three dimensional 0(2) model. Renormalization 
group arguments predict a mild non-analyticity in the specific heat, c v ~ t~ a with 
t — (T — T c )/T c and a = —0.0151(3) |43j. This prediction agrees reasonably well with 
micro gravity experiments which find a = —0.01285(4) |44j . 

The excitation spectrum of superfluid 4 He shows important differences as compared 
to the spectrum of a dilute Bose condensed gas. As expected, at small momenta the 
excitations are phonons with a linear dispersion relation e(p) = c s p, where the speed 
of sound at normal pressure is c s = 238 m/sec. At larger momenta the dispersion 
relation has a second minimum, called the roton minimum. The dispersion relation in 
the vicinity of the minimum is 

«p)=a + ^, ^ («» 

where A/ks = 8.7 K, m* = 0.14 m andpo/h = 1.9 A 1 . The roton plays a significant role 
in determining the specific heat and transport properties near the critical temperature. 
The existence of the roton is closely related with strong short range correlations in liquid 
helium. These correlations can be quantified in terms of the static structure factor S(q), 
which is the Fourier transform of the density correlation function 

S(q) = -j d*xe~^ [(p(O)p(x)) - (p(0)) 2 ] . (17) 

The static structure factor vanishes linearly in q for small momenta, and approaches 
a constant value for large q. S(q) has a sharp maximum at intermediate values of 
q, which reflects the presence of correlations at the scale of the typical inter-atomic 
distance. Feynman proposed a variational wave function for excitations in liquid helium 
which gives e(q) = q 2 /(2mS(q)) [45]. This relation can also be derived from an effective 
hydrodynamic Hamiltonian, see [25]. Feynman's result reproduces both the phonon 
dispersion relation at low momentum, and the roton minimum at larger momentum. 
In order to study the kinetics of liquid helium one has to understand the scattering of 
phonons and rotons. The phonon-phonon and phonon-roton interaction is determined 
by the equation of state and by constraints from Galilean invariance and U(l) symmetry 
|25j . We will discuss these constraints in more detail in the next section, in connection 
with superfluid Fermi gases. 

2.3. Fermi liquids: The dilute Fermi gas at unitarity 

In this section we consider non-relativistic Fermi liquids. Fermionic systems are 
interesting because it is possible to make strongly correlated liquids with only zero 
range interactions. The Fermi liquid is described by the same lagrangian as in equ. (jE]) 

£ = tf (id + ^ V - y (V^) 2 , (18) 
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were ip is now a two-component fermion field with mass m. The coupling constants 
Cq is related to the scattering length by the same relation as in the bosonic case, 
C = (ATih 2 a)/m. 

Over the last ten years there has been truly remarkable progress in the study of 
cold, dilute gases of fermionic atoms in which the scattering length a of the atoms 
can be controlled experimentally. These systems can be realized in the laboratory 
using Feshbach resonances, see jl6] for a review. A small negative scattering length 
corresponds to a weak attractive interaction between the atoms. This regime is known as 
the BCS (Bardeen-Cooper-Schrieffer) limit. As the strength of the attractive interaction 
increases the scattering length becomes larger. It diverges at the point where a two- 
body bound state is formed. The point a = oo is called the unitarity limit, because 
the scattering cross section saturates the s-wave unitarity bound a = 4ir/k 2 . On the 
other side of the resonance the scattering length is positive. For large positive values 
of a the two-body binding energy is related to the scattering length by B = h 2 /(ma 2 ). 
The regime in which the binding energy becomes large is called the BEC (Bose-Einstein 
condensation) limit. 

We now consider properties of the the many-body system as a function of the s- 
wave scattering lengths. In the high temperature limit the equation of state is again 
that of an ideal gas, and the leading correction is described by the virial expansion. For 
small a the second virial coefficient is 

' =A 3 + -aA 2 . (19) 



8V2 2 

In the limit a — > oo the interaction term goes to — X s / (2^/2). The Fermi gas becomes 
degenerate as nA 3 ~ 1. In the limit in which the scattering length is large the Fermi gas 
becomes strongly interacting at the same temperature at which quantum effects become 
important. 

At low temperature and in the BCS limit, a < and n 1 / 3 ^! < 1, the Fermi gas can 
be described as a Landau Fermi liquid. The excitations are weakly interacting particles 
and holes which carry the quantum numbers of the elementary fermions. At very low 
temperature the particle-particle interaction near the Fermi surface becomes large, and 
the Fermi liquid undergoes a phase transition to a BCS superfluid. The transition 
temperature is [UJ 

8e^'E F / i \ 

Tc= (4iF^% eXP (~2^Hj' (20) 
where 7 is the Euler constant. The Fermi momentum kp is defined by the relation 

h.3 

<"> 

and Ep = k F /(2m) is the Fermi energy. This relation defines a "Fermi momentum" even 
in the case that no sharp Fermi surface exists. Note that Tp = Ep is the degeneracy 
temperature (we have set ks = 1). Also note that rz. 1 / 3 1 a | < 1 implies T c <^Tp. 

In the Bose-Einstein limit the fermions form tightly bound molecules. The residual 
interaction between the molecules is repulsive, and the many-body system behaves as 
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a weakly non-ideal Bose gas. The Bose gas condenses at the critical temperature given 
in equ. ([TO]) . Using the fact that the mass of molecules is 2m, and that their density is 
n/2, we get 

T c = 0.21E F . (22) 

Variational calculations suggest that at zero temperature the evolution from weak to 
strong coupling is smooth [J5] . The system is a pair condensate for all values of the 
coupling, but the size of the pairs evolves from being much smaller than the inter-particle 
spacing in the BEC limit to being much larger in the BCS limit. This idea is confirmed 
by quantum Monte Carlo calculations [19] and experimental observations |50j . 

Of particular interest is the crossover ( "unitarity" ) regime where a — > oo. The 
Fermi gas at unitarity possesses a number of interesting properties. First of all, the 
system is scale invariant [5TJ[52]. This implies, for example, that all energy scales in the 
many body system, like the critical temperature, the gap, and the chemical potential, 
are proportional to the Fermi energy 

T c = aE F , A = PE F , fi = £E F . (23) 

Similarly, all length scales are given by numerical constants times the inverse Fermi 
momentum. The values of the universal constants a, /?,£,.. . can be determined using 
Quantum Monte Carlo (QMC) calculations, or from experiments on harmonically 
trapped Fermions. QMC calculations performed by Burovski et al. give T c = 0. 152(7) E F 
[53], and Carlson et al. obtained A = 0.50(3)£ F [54] and /z = 0AA(1)E F [55]. A 
summary of experimental results was recently given by Luo and Thomas [56] . 

Second, the unitarity regime is the most strongly correlated simple many body 
system. The crossover regime is continuously connected to both the non-interacting 
Fermi gas and the non-interacting Bose gas, but neither limit provides a quantitatively 
accurate description. Very important for the purpose of this review is the observation 
hydrodynamic behavior and low viscosity in very dilute Fermi gases in the unitarity 
limit. 

In order to study the kinetic description of a dilute Fermi gas at unitarity we 
have to determine the nature of the quasi-particles and their interaction. In the high 
temperature limit the excitations are elementary fermions, even in the limit a — > oo. 
This is related to the fact that the average cross section is of order A 2 , where A ~ T -1 / 2 
is the thermal wave length. In the low temperature superfluid phase the dominant 
excitation is the phonon. The dispersion relation is 

e p = c s p, c s = J| v F , (24) 

where v F = k F /m is the Fermi velocity, and £ is the universal parameter defined in 
equ. (1231) . Corrections to equ. ff24l) are of the order p 2 /(m/i) [57], and become large 
when p ~ hp. The dispersion relation for momenta near k F is not well constrained. The 
static structure factor has been measured in quantum Monte Carlo simulations, and it 
does not show a liquid-like peak [58] . This suggest that the dispersion relation does not 
have a roton minimum. 
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The three and four-phonon interaction is completely fixed by the equation of state 
and symmetry constraints. These constraints are most easily derived from an effective 
lagrangian for the phonon field. The phonon field is defined as the phase of the order 
parameter 

(W) = \(W)\e 2Ufi . (25) 

We now construct the most general lagrangian for the field p which is consistent with 
Galilei invariance and U(l) symmetry. A U(l) transformation changes the phase of 
the wave function and acts as a shift on phonon field, p — *> ip + a. Invariance under 
the U(l) symmetry requires that the lagrangian only contains derivatives of p. The 
phase symmetry can be extended to time-dependent transformations xp —>■ exp(ia(t))tp 
if the chemical potential transform as fi — > fi + Tid^a. This is a symmetry of the 
effective lagrangian if the chemical potential always appears in the combination fi+hdop. 
Finally, under a Galilei transformation with velocity v the phonon transforms as 
y?(x, t) — > </?(x — v£) — mv • x + 0(v 2 ). This implies that time derivatives of p have to 
be accompanied by spatial derivatives of (p. At leading order in derivatives of p we can 
incorporate these symmetries by constructing a lagrangian that only depends on the 
variable 

X = ll - m „-^. (26) 

The functional form of the lagrangian C(X) is fixed by the observation that for a 
constant phonon field the lagrangian reduces to a function of \i. Since differentiating the 
lagrangian with respect to the chemical potential gives the density this function must 
be the pressure -P(/i). We conclude that 

where we have used the fact that, up to a numerical factor, the pressure of the interacting 
system is equal to that of a free gas. We have also used that this factor can be related 
to the ratio £ = iijEp. Phonons are low energy excitations and we can expand equ. ( 1271) 
in powers of dop and Vip. We find 

£ = kdo<P? - k ( W) 2 - « [(<9 O 0) 3 - 9c^ o (V0) 2 l (28) 



2 V ' 2 

3 2 
a 

2 



(9 o 0) 4 + 18^(a o 0r(V0) 2 -27 C 4 



+ 



where we have rescaled the field p = const x to make it canonically normalized. We 
have also defined a = 7rc^ 2 £ 3 / 4 /(3 1 / 4 8// 2 ). We observe that the three and four phonon 
vertices are completely fixed by the speed of sound c s . This implies that there are 
no free parameters that enter into the kinetic theory of phonons. We also note that 
equ. 02 71) generates phonons self interactions to arbitrary order in the phonon field, but 
to leading order in the number of derivatives. Terms involving higher derivatives of p 
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were constructed in [52] • These terms involve non-trivial constraints from not just scale 
invariance, but from the full conformal symmetry of the Fermi gas at unitarity. 

About Units: Up to this point, we have explicitly displayed factors of h, c and ks- From 
now on we will work in natural units and set h = = c = 1. 



2-4- Gauge theories: QCD 

Quantumchromodynamics (QCD) is governed by the lagrangian 

£ = —Gl v G% + £ ^ - m f )1>f , (29) 

where ipf is a Dirac fermion with flavor index / and my is the quark mass. We have 
suppressed the color (A = 1, . . . , N c ) and spinor (a = 1, . . . , 4) indices of the fermion 
fields. The covariant derivative acting on the quark fields is 

z^ = 7 M (^ + ^y^, (30) 

where A a ^ is a gauge potential and A a (a = 1, . . . , N% — 1) are the Gell-Mann matrices. 
The field strength tensor is defined by 

G% = d,Al - d v Al + gr hc A\Al , (31) 

where f abc are the SU(N C ) structure constants, and g is a coupling constant. In the 
standard model N c = 3 and Nj = 6, but three out of the six flavors are too heavy too 
play much of a role in the dynamics of QCD, and we shall mostly concentrate on Nf = 3 
flavors. The total quark density 

P q = H 1>tyf ( 32 ) 
/ 

is conserved and we can introduce a chemical potential fi coupled to p q . The phase 
structure and transport properties of QCD at finite \i are an interesting subject |59j, 
but in this review we will concentrate on QCD at non-zero temperature and zero or 
very small small chemical potential. It is interesting to note that at low quark density 
the relevant degrees of freedom are protons and neutrons. In the low energy limit 
the interaction between neutrons and protons is governed by an effective lagrangian of 
the type given in equ. (JTSj). The scattering length is a function of the quark masses, 
and it is theoretically possible to tune the light quark masses to a point where the 
neutron-neutron scattering length diverges. The real world is close to this point, as the 
experimental value of the scattering is a nn ~ — 17 fm is much larger than typical QCD 
scales. This implies that there is a point in the QCD phase diagram where the long 
distance physics is equivalent to that of a dilute atomic Fermi gas at unitarity. 

For many purposes we can consider the first three flavors (up, down, and strange) 
to be approximately massless. In this limit the QCD lagrangian contains a single 
dimensionless parameter, the coupling constant g. If quantum effects are taken into 
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account the coupling becomes scale dependent. At leading order the running coupling 
constant is 

9 2 {q) = 7 j S , 6o = ^e-^/. (33) 

6 log(g 2 /A^ CT ) 3 3 

This result implies that as a quantum theory, QCD is not characterized by a 
dimensionless coupling, but by a dimensionful scale, the QCD scale parameter Aq C d- 
This phenomenon is called dimensional transmutation. We also observe that the 
coupling decreases with increasing momentum. This is the phenomenon of asymptotic 
freedom. 

At high temperature the dominant momenta are on the order of T, and for 
T ^> Aqcd asymptotic freedom implies that bulk thermodynamics is governed by weak 
coupling. The weak coupling expansion of the equation of state is 

P = T 4 {c + c 2 g 2 + c 3 g 3 + (c' 4 \og(g) + c 4 )g 4 + ...}, (34) 

where the first term is the Stefan-Boltzmann law and 

Co = ^(2(A 6 2 -l)+4A c A / -) (35) 

depends on the number of degrees of freedom (2 (A 2 — 1) gluons and 4A C A/ quarks). 
We note that in a theory of massless particles the equation of state is always sensitive 
to quantum statistics, even if the temperature is high. The first correction is [60] 



/V 2 — 1 / 5 \ 

" = ~iM-( N ' + i N ')- (36) 

The perturbative expansion in equ. (l34"j) is evaluated with g taken to be the running 
coupling constant evaluated at a scale q ~ T. The precise scale is not uniquely 
determined - changing the scale corresponds to reshuffling higher order corrections in 
the perturbative expansion. The scale is usually chosen to improve the apparent rate of 
convergence. This criterion gives a value close to 2nT . 

We note that the perturbative expansion is not a power series in the fine structure 
constant a s = g 2 / (4?r). The expansion contains square roots and logarithms of a s . Non- 
analytic terms in the expansion are related to infrared sensitive diagrams. For example, 
the g 3 term in equ. (134"I) is due to ring diagrams (also called the plasmon term). Ring 
diagrams are one-loop gluon diagrams in which the leading order gluon self energy has 
been summed to all orders. We also note that the weak coupling expansion cannot be 
extended to arbitrarily high powers in g. At 0(g 6 ) one encounters infrared divergent 
diagrams which can only be summed non-perturbatively, by computing the partition 
function of three-dimensional QCD at zero temperature. 

In order to analyze the relevant scales in high temperature QCD in more detail we 
consider the current- current interaction 

M = jfil^jl , (37) 

where is a color current and 11^ is the gluon polarization function. The tensor 
structure of the gluon polarization function can be decomposed into a transverse and a 
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longitudinal part 

= Tf{q)Pl + U L (q)P^ (38) 



We will consider the polarization function in the limit of weak coupling [g < 1), and for 
uj ^ q ^ T, where uj, q are the energy and momentum transfer. We find 

sab 

n afei (g) = ^— — , (39) 

sab 

n afeT (<?) = 2 % 2 5 , (40) 



where 

m 2 D = ^ 2 T 2 (l + ^ j (41) 

is called the Debye mass. The longitudinal term governs the color-Coulomb interaction 
between static charges. We observe that the Coulomb interaction is screened at 
distances r ~ m^ 1 ~ l/(gT). In perturbation theory the static magnetic interaction 
is unscreened [16], but non-static magnetic interactions are dynamically screened at 
a distance r ~ {m 2 D uj)~ l l z . This phenomenon, known as Landau damping, is due to 
the coupling of gluons to particle-hole (or particle-anti-particle) pairs, and also play 
a role in electro-magnetic plasmas. Unlike classical plasmas the QCD plasma has a 
non-perturbative static magnetic screening mass tum ~ g 2 T. This is the scale that 
determines the non-perturbative g 6 term in the pressure. Modes below the magnetic 
screening scale contribute 

P ~ T / d 3 k ~ g 6 T 4 . (42) 



The gluon polarization tensor also determines the propagation of gluonic modes. For 
this purpose we need the full energy and momentum dependence of U T ' L , see [61] . For 
momenta q ^> gT there are two transverse modes with dispersion relation uj ~ q. For 
momenta q < gT there are two transverse and one longitudinal mode. The longitudinal 
mode is sometimes called the plasmon. The energy of all three modes approaches 
uj = Up = mc/v3 as q — > 0. The quantity uj p is known as the plasma frequency. The 
gluon (and plasmon) decay constant in the limit q — > is [62] 



7 = 6.64^^ . (43) 
' 24vr v ; 

An important issue is how small the coupling has to be in order for the perturbative 

estimates to be applicable. The convergence properties of the weak coupling expansion 

for the pressure are extremely poor. The series shows no signs of converging unless the 

coupling is taken to be much smaller than one, j < 1, corresponding to completely 

unrealistic temperatures on the order of 1 TeV. The problem is mostly due to the 

non-analytic terms in the expansion, and convergence can be improved significantly 
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by considering resummation schemes or self-consistent quasi-particle expansions 
Convergence can also be improved by using a hierarchy of effective field theories for the 
hard (p ~ T), electric (p ~ gT), and magnetic {p ~ g 2 T) sectors of the QCD plasma 
[53] . Ordinary perturbation theory corresponds to treating the hard and the electric 
sector perturbatively, but convergence can be improved by treating both the electric 
and the magnetic sector non-perturbatively [65] . 

Despite these advances accurate results at temperatures that can be reached in 
heavy ion collisions at RHIC have to rely on numerical simulations of the QCD partition 
function on a space-time lattice, see [66] for a review. Lattice simulations with realistic 
quark masses find a phase transition at the critical temperature T c = 192(8) MeV [67J. 
The transition is a rapid (but smooth) crossover from a low temperature phase that 
exhibits chiral symmetry breaking and confinement to a chirally restored and deconfined 



high temperature phasej. The energy density reaches about 85% of the ideal gas value 
at T ~ 2T C and then evolves very slowly towards the non- interacting limit. 

Below the critical temperature the degrees of freedom are hadrons. The lightest 
hadrons are pions, which are the Goldstone bosons associated with the spontaneous 
breaking of the chiral symmetry of the QCD lagrangian. We can view pions as a 
spin-isospin sound wave that propagates in the QCD vacuum. Because quarks are 
not massless the chiral symmetry is not exact, and pions have non- zero masses. The 
masses of the charged and neutral pions are m n ± = 139 MeV and m n o = 135 MeV. The 
lightest particle which is not a type of sound wave is the rho meson, with a mass of 
770 MeV. Chiral symmetry constrains the pion scattering amplitudes. As in the case of 
phonons, these constraints are obtained most easily from the low energy effective chiral 
lagrangian. At leading order we have 



f 2 

£ = ^Tr 
4 



c^Ec^E+l + \BTr(M^) + h.c] + . . . , (44) 



where £ = exp(i(f> a X a / f v ) (a = 1, ... ,8) is the chiral field, f n = 93 MeV is the pion 
decay constant, B is proportional to the quark condensate, and M = diag(m u , m^, m s ) 
is the mass matrix. We note that f n can be viewed as the stiffness of the QCD vacuum, 
?m f)P 
om q 

where P vac ~ 0.5 GeV/fm 3 is the vacuum pressure, and m q = {m u + m d )/2. This result 
follows from the Gell-Mann-Oakes-Renner relation m 2 f% = (m u + m^) (iptp) together 
with (dP vac )/(dm) = (ipijj)- An expansion of S in powers of the field <j) a determines the 
interaction between pions. Restricting ourselves to the SU{2) flavor sector (pions only) 
we get 

C = \{d,rf - \m 2 Mf + [( W 1 ) 2 - m\d,<p b f] + . . . , (46) 

* This issue is not completely settled. Aoki et al. find distinct crossover transitions at significantly lower 
temperatures, T x — 151 MeV for chiral symmetry restoration, and Ta.ec = 175 MeV for deconfincmcnt 



Nearly Perfect Fluidity 



18 



s / s ° Pure S lue Q CD S/S a N = 4 super- Yang-Mills 

7 6 5 MiiMS=2-*T 




- weak-coupling to order A 3,/2 
-strong-coupling to order A -3 / 2 



T/T c 




Figure 1. Entropy density in units of the Stefan-Boltzmann value for pure gauge QCD 
and Af — 4 supersymmctric QCD. The left panel shows the entropy density of pure 
gauge QCD as a function of T/T c . The grey band is the lattice result. The solid lines 
show a resummed QCD calculation 63J. The different lines correspond to different 
choices for a non-perturbative parameter ca- The dashed lines mark an error band 
determined by variations in the QCD renormalization scale. The right panel shows the 
entropy density of SUSY QCD as a function of the 't Hooft coupling A. The curves 
arc labeled as in the left panel. 



where cj) a (a = 1, 2, 3) is the pion field. This result is clearly analogous to the phonon 
interaction in equ. (I28I) . There are, however, some minor differences. Because of parity 
and isospin symmetry there are no vertices with an odd number of pions. We also note 
that the leading four-pion interaction has two derivatives, while the four-phonon term 
involves four derivatives. 

2.5. Gauge theories: Superconformal QCD 

QCD is a complicated theory, and a significant amount of effort has been devoted to the 
study of generalizations of QCD that possess a larger amount of symmetry, in particular 
super symmetry. Supersymmetry is a symmetry that relates bosonic and fermionic fields. 
The simplest supersymmetric cousin of QCD is SUSY gluodynamics, a theory of gluons 
and massless fermions in the adjoint representation of the color group called gluinos. 
Theories with more supersymmetry involve extra fermions and colored scalar fields. The 
most supersymmetric extension of QCD is a theory with four supersymmetries, called 
M = 4 SUSY QCD. Theories with even more supersymmetry contain fields with spin 
3/2 and 2, and therefore involve gravitational interactions. These theories are known 
as supergravity. 

The lagrangian of Af = 4 SUSY QCD is 

C = -\Gl v G% - fXy^D^ + D^D^% + £ xx<t> + , (47) 

where G a is the usual field strength tensor, A" is the gluino field, and (f>^ is a colored 
Higgs field. The gluino is a two-component (Weyl) fermion in the adjoint representation 
of the color group. The index i (i = 1,...,4) transforms in the fundamental 
representation of a global SU(4)r "R-symmetry" . This symmetry interchanges the 
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bosons and fermions that are related by the four supersymmetries, and is analogous to 
the flavor symmetry of QCD. The Higgs is a scalar field in the adjoint representation of 
color, and in an anti-symmetric tensor (six dimensional) of SU(4)r. Note that the total 
number of fermionic fields, 8(iV| — 1), is indeed equal to the number of bosonic fields. 
We have not explicitly displayed the Yukawa couplings £aa</> and Higgs self couplings 
£^4, see [69]. Both interaction terms only involve the dimensionless gauge coupling g. 

M = 4 SUSY QCD has a vanishing beta function and is believed to be a 
conformal field theory (CFT). As a consequence there is no dimensional transmutation, 
no confinement or spontaneous symmetry breaking, and no phase transition. The theory 
is an Coulomb phase for all values of the coupling g and the temperature T. However, 
if g is not small then there is no obvious way to compute thermodynamic or transport 
properties of the plasma. 

An interesting new approach is provided by the duality between strongly coupled 
large N c gauge theory and weakly coupled string theory on AdS 5 x S 5 discovered by 
Maldacena [20] . We will have more to say about this approach in Section HI For now 
we observe that the correspondence can be extended to finite temperature. In this case 
the relevant configurations is an AdSs x S5 black hole. The temperature of the gauge 
theory is given by the Hawking temperature of the black hole, and the entropy is given 
by the Hawking-Beckenstein formula S = A/ (AG), where A is the surface area of the 
event horizon and G is Newton's constant. 

The AdS/CFT correspondence makes predictions for the thermodynamics of the 
gauge theory in the limit of a large number of colors, N c — > 00. The perturbative 
expansion of a SU(N C ) Yang-Mills theory involves the 't Hooft coupling A = g 2 N c . 
In the weak coupling limit we take N c — > 00 with A = const and A < 1. Using 
the AdS/CFT correspondence we can also study the strong coupling limit N c — *■ 00 
with A ^> 1. M = 4 SUSY QCD is a conformal theory and scale invariance implies 
that e = 3P as well as s = 4P/T, where e is the energy density and s is the entropy 
density. By dimensional analysis the entropy density of the interacting system must be 
proportional to the entropy density so of the free system. The weak and strong coupling 
expansions for s/s are [701 [7TI [72] 



This result has a number of remarkable features. First we observe that the entropy 
density at infinite coupling only differs by a factor 3/4 from the result in the free theory. 
We also note that the first non-trivial corrections in the strong and weak coupling 
limit are consistent with the idea that the evolution from weak to strong coupling is 
smooth. Equation (1481) was compared with resummed perturbation theory and Pade 
approximants in [73], see Fig. [TJ The authors argue that at the "QCD-like" point 
s/sq = 0.85 neither the strong nor the weak coupling expansion are quantitatively 
reliable, but that resummed perturbation theory is useful in this regime. 




A > 1, 
A < 1. 



(48) 
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3. Transport theory 

In this section we summarize theoretical approaches to transport phenomena in strongly 
coupled quantum fluids. The most general of these approaches is hydrodynamics. 
Hydrodynamics is based on the observation that correlation functions at low energy 
and small momentum are governed by the evolution of conserved charges. Conservation 
laws imply that the densities of conserved charges cannot relax locally, but have to 
propagate or diffuse out to large distance. This corresponds to hydrodynamic excitations 
with dispersion laws of the form u ~ q (sound) or uj ~ iq 2 (diffusion). 

Hydrodynamics can be developed as an expansion in derivatives of the fluid velocity 
and the thermodynamic variables. The leading order theory, called ideal hydrodynamics, 
only depends on the equation of state, and is exactly time-reversible. The next order 
theory, (first order) viscous hydrodynamics, involves a new set of parameters called 
transport coefficients, and describes dissipative, time-irreversible phenomena. The 
validity of hydrodynamics is controlled by the relative size of the next-to-leading order 
terms. If dissipation is dominated by shear viscosity§ then the expansion parameter is 
I /Re, where Re is the Reynolds number defined in equ. (j4j). 

The values of the transport coefficients can be extracted from experiment, or 
computed from an underlying field theory. The connection between transport coefficients 
and correlation functions in a (quantum) field theory is provided by linear response 
theory. Using linear response theory one can relate transport coefficients to the zero 
energy and zero momentum limit of a retarded correlation function. These relations are 
known as Kubo formulas, see Sect. 13.41 

Calculations based on the Kubo formula are difficult, in particular if the interaction 
is not weak. The situation simplifies if the system allows a microscopic description in 
terms of quasi-particles. In that case we can use an intermediate effective theory, known 
as kinetic theory, to relate the microscopic lagrangian to the hydrodynamic description. 
Kinetic theory also provides a more microscopic criterion for the applicability of 
hydrodynamics. Using the kinetic estimate for the shear viscosity in equ. (j2j) we get 
1/Re ~ (v/c s )(l m fp/L) where c s is the speed of sound and the ratio 

Kn = 1 ^ (49) 

LJ 

is called the Knudsen number. Hydrodynamics is valid if the mean free path is much 
smaller than the characteristic size, and the Knudsen number is small. 

The calculation of transport coefficients in kinetic theory is reviewed in Sect. 13.51 
If the interaction between quasi-particles is strong then the kinetic description breaks 
down. A new approach to extracting transport properties from a strongly coupled field 
theory is the holographic method which we will discuss in the next section. Using 
holography the calculation of the retarded correlator can be reduced to a classical 
computation in a suitable dual theory. 

§ If dissipation is dominated by heat transport, then the expansion parameter is l/(Re ■ Pr), where 
Pr is the Prandtl number defined in equ. (I140j) . 
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3.1. Hydrodynamics 

3.1.1. Non-relativistic fluids The hydrodynamics of a one-component non-relativistic 
fluid is governed by the conservation laws of energy, mass (particle number), and 
momentum, 

§ + V-r =0, (50) 
^ + V-g =0, (51) 

^ + V,n 4 , = . (52) 

Here, e is the energy density, p is the mass density, g is the momentum density, and H^- is 
the stress tensor. The relations between the conserved currents and the hydrodynamic 
variables are called constitutive relations. These relations can be determined order 
by order in an expansion in derivatives of the flow velocity and the thermodynamic 
variables. The leading order result is called "ideal hydrodynamics". At this order the 
constitutive relations are completely fixed by Galilean invariance, rotational invariance, 
and conservation of entropy. The result is 

j e =v(e + P), (53) 
g = pv , (54) 
Uij = P5ij + pViVj , (55) 

where e = eo + \pv 2 and eo is the energy density in the rest frame of the fluid. There are 
six hydrodynamic variables, v, p, e and P, which are determined by the five conservation 
laws fl50H52p . In order for the equations to close we need to supply an equation of state 
P = P(e,p). Since hydrodynamic variables evolve slowly, the equation of state is the 
one in thermal equilibrium. 

In ideal hydrodynamics the equations of continuity and momentum conservation 

are 

f + V-(pv) =0, (56) 
<9v 1 

_ + (v V)v= --VP. (57) 

The equation of momentum conservation is known as the Euler equation. In the 
case of ideal hydrodynamics the equation of energy conservation can be rewritten as 
conservation of entropy, 

^ + V-(v S ) = 0. (58) 

At next order in the derivative expansion dissipative terms appear. The size of these 
terms is controlled by new parameters called transport coefficients. The relation g = pv 
is not modified (it follows from Galilean invariance), but two new coefficients appear in 
the stress tensor. We can write = P5ij + pvtVj + d~Uij with 

511^ = -77 (v^ - VjVi - • v) - (8 i:j (V • v) . (59) 
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Here, r\ is the shear viscosity and £ is the bulk viscosity. The correction to the energy 
current has the form j\ = t>j(e + P) + VjdHij + Qi with 

Q = -kVT, (60) 

where T is the temperature and k is the thermal conductivity. The second law of 
thermodynamics implies that rj, (, k > 0. The equation of momentum conservation with 
the viscous stresses (159]) included is known as the Navier-Stokes equation. 

The linearized hydrodynamic equations describe the propagation of sound and 
diffusive modes. In the case of a non-relativistic fluid there is a pair of sound modes 
that couple to the pressure/density and the longitudinal velocity, a pair of diffusive 
shear modes that couple to the transverse velocity, and a diffusive heat mode. The 
longitudinal and transverse components of the velocity are defined by v = v L + v T 
with Vxv t = and V • v T = 0. The hydrodynamic modes govern the hydrodynamic 
correlation functions. The transverse velocity correlation function is defined by 

S%(oj,k) = (SvfSvj)^ = J d 3 xdt e l(ujt -^(5vf(x,t)5vj(0,0)} , (61) 

where 5vf(x,t) = t>f(x, t) — (ff(x, t)) is a fluctuation of the velocity. Linearized 
hydrodynamics gives 

^(^ k ) = y(^-^')^^ki' ( 62 ) 

where v = rj/p is the kinetic viscosity. The dependence of the correlation function on 
uj and k is determined by the laws of hydrodynamics, equ. (I50H52I) . and the overall 
normalization is fixed by the thermodynamic relation 

T 

(5fj(x, t)5vj(-x, t)) = —5ij5(x. — x') . (63) 

We observe that the transverse velocity correlation function has a diffusive pole, where 
the diffusion constant is given by the kinematic viscosity. The entropy correlation 
function has a diffusive pole governed by the thermal diffusion constant \ = k/(c p p), 
where c p is the specific heat at constant pressure. The correlation function is 

S"(u, k) = (8s5s)^ = 2 X . (64) 

The pressure correlation function contains the sound pole and is given by 

S pp {u, k = (6pSp) wM = 4pTc s — , 65 

where c s = [{dP)/(dp)\ s \ 1 ^ 2 is the speed of sound and 7 = 7^ + 7t is the coefficient of 
sound absorption (the inverse sound attenuation length). The contributions to 7 from 
viscosity and thermal conductivity are given by 

k ' fr-M 4 \ ( dT \ 

7 ^ = 27^l c + 3^' lT = ^r x \dp) s - (66) 

These results illustrates the criterion for the validity of hydrodynamics given above. 
Hydrodynamics is based on a small momentum expansion. Applied to equ. fl6"5"l) this 
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implies that uo ~ c s k ^> k 2 i]/p. Taking the characteristic size to be L ~ 1/k this 
is equivalent to r]/(c s pL) < 1, A more microscopic criterion follows from the kinetic 
estimate of the shear viscosity given in equ. (j2J): Hydrodynamics describes sound waves 
with a wave length that is large compared to the mean free path: L ^> l m f p . 

3.1.2. Superfluid hydrodynamics Superfluidity is characterized by the spontaneous 
breakdown of the U(l) symmetry associated with the conserved particle number. By 
Goldstone's theorem the spontaneous breaking of a continuous symmetry leads to the 
appearance of a gapless mode. This mode has to be included in the hydrodynamic 
description of the system. We introduced the Goldstone boson field <p in equ. (|25|) . 
The quantity v s = V<p/ m can be interpreted as the superfluid velocity. Since v s is the 
gradient of a phase the superfluid velocity is irrotational, V x v s = 0. 

We have to generalize the constitutive equations to include both the normal fluid 
velocity v n and the superfluid velocity v s . In the ideal fluid case (no dissipation) the 
result is completely fixed by Galilean invariance and thermodynamic relations. The 
constitutive equations are 

g = PnV„ + p s v s , (67) 
Uij = Pdij + p n v n ,iV n j + p s v s ,iVs d , (68) 

j e = psTv n + ffi + -vfj (p n v n + p s v s ) + p„v„v„ ■ (v n - v s ) , (69) 

where p n and p s are the normal and superfluid density of the system. The total density 
P = Pn + Ps is the sum of the normal and superfluid contributions. The ratio p s /p is a 
function of the temperature, the chemical potential, and the relative velocity |v n — v s |. 
This function, like the equation of state P(p,T,\v n — v g |), depends on microscopic 
details. The conservation laws are given by equ. (I50H52I) . These equations have to be 
supplemented by an equation of motion for the superfluid velocity. Landau showed that 
Euler's equation for the superfluid velocity is given by [74] 

^ + (v s .V)v s = -V/i. (70) 

Because v s is irrotational we can rewrite the convective derivative on the LHS of 
equ. (17D1 as a total derivative, v s ■ Vv s = |V(v s 2 ). 

As in the case of a normal fluid we may consider dissipative corrections to the 
constitutive equations. The form of these terms is constrained by rotational and Galilean 
invariance, and by the second law of thermodynamics. The viscous corrections to the 
energy momentum tensor are 

5Ylij = - 7] (ViV nJ + VjV^i - g^ijV • V n ^j (71) 



- <%(CiV ■ (p s (v s - v n )) + C2 (V • v r , 

We observe that viscous shear stresses only arise from the normal component of the 
flow. In addition to the normal bulk viscosity term proportional to (2 there is a second 
contribution that involves the relative motion of the normal and superfluid components. 
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Two additional bulk viscosities appear in the dissipative correction to the RHS of 
equ. (1701) . We replace V/x by V(/i + H) with 

H = -CsV • (p s (v, - v n )) - C 4 V • v n . (72) 

Onsager's symmetry principle requires that £ 4 = The dissipative correction to the 
energy current is bj\ = v n jSHij + p s (v s ,i — v Ht i)H + Qi where Qi = — kVjT as in the case 
of a normal fluid. 

Superfluid hydrodynamics contains two velocity fields, the normal flow velocity v n 
and the superfluid (irrotational) flow velocity v s . This extra degree of freedom leads to 
an additional sound mode called second sound. The velocity of second sound depends 
strongly on temperature and vanishes at the critical temperature where p s /p — > 0. 
If thermal expansion can be neglected second sound is an oscillatory motion of the 
superfluid against the normal fluid which does not lead to any mass transport and can 
be viewed as a pure entropy wave. 



3.1.3. Relativistic fluids In a relativistic fluid the equations of energy and momentum 
conservation can be written as a single equation 

8^ = 0, (73) 

where T^ v is the energy momentum tensor. In ideal fluid dynamics the form of is 
completely fixed by Lorentz invariance, 

= (e + P)u^u v + P^ u , (74) 

where is the fluid velocity (u 2 = —1) and vf v = diag(— 1, 1, 1, 1) is the metric tensor. 
In a relativistic theory there need not be a conserved particle number. If a conserved 
particle number, for example baryon number, exists then there is a second hydrodynamic 
equation that express particle number conservation 

<9 M (n<) = 0, (75) 

where n is the particle density. As in the non-relativistic case the hydrodynamic 
equations have to be supplemented by an equation of state P = P(e) or P = P(e,n). 
The four equations given in equ. (1731 can be split into two sets using the longitudinal 
and transverse projectors 

A[L = ~ u i^ u u, = r}^ + u^Uy . (76) 

With the help of the thermodynamic relations de = Tds and e + P = sT the longitudinal 
equation is equivalent to entropy conservation 

d M («O = 0, (77) 

and the transverse equation is the relativistic Euler equation 

Du^-j^ViP, ( 78 ) 

where D = u ■ d and V^; = A^ u d u . Comparison with equ. (1571) shows that the inertia of 
a relativistic fluid is governed by e + P. 
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The form of the dissipative terms depends on the precise definition of the fluid 
velocity. A useful choice is to define by the requirement that in the local rest frame 
T 00 = e and T° l = 0. This definition is called the Landau frame [74J. In this frame 
the dissipative correction to the energy momentum tensor in the rest frame has the 
same form as in the non-relativistic case, see equ. (1591) . We write the stress tensor as 

= T<f + 5^T^ + 5 {2) T^ + . . ., where Tg v is the stress tensor of the ideal fluid given 
m equ. flUD, 6 W T^ is the first order viscous correction, etc. A covariant expression for 
SWT^ is 

S (i) T ^ = _ va ^ _ ^ A ^ d . u ( 79 ) 
where we have defined 



2 



A^A^ d a u p + d p u a - -r, aP d ■ u . (80) 



~ ~ ""3" 
The dissipative correction to the conserved particle current is j M = nu^ + 5j M with 

^ = -«(^)'<£). (81) 
where k is the thermal conductivity and /i is the chemical potential associated with the 
conserved density n. Alternatively, one can define the velocity via the conserved particle 
current (Eckart frame). In that case there is no dissipative contribution to j M and the 
thermal conductivity appears in stress tensor. 

The hydrodynamic equations determine the propagation of sound and diffusive 
modes. We consider the case without a conserved particle number. In this case all 
the modes can be found by considering correlation functions of the energy-momentum 
current g l = T°\ The longitudinal and transverse correlation functions are 

I> 2 k 2 

{u 2 - c^k 2 ) 2 + (r s uk 2 y 



99 



As in the non-relativistic fluid we find a pair of sound waves, and a pair of diffusive 
shear modes. The sound attenuation length is given by 

Ts = , (84) 

and the analog of the kinematic viscosity is the ratio r]/(sT). 

An new issue that arises in viscous relativistic hydrodynamics is the apparent lack 
of causality of the equations of motion. The problem can be seen by inspecting the 
linearized equation for the diffusive shear mode. The equation is first order in time, but 
second order in spatial gradients. As a result discontinuities in the initial conditions 
can propagate with infinite speed. This is not really a problem of the hydrodynamic 
description - the relevant modes are outside the domain of validity of hydrodynamics 
- but the acausal modes cause difficulties in numerical implementations. To overcome 
these difficulties one can include second order gradient corrections in the stress tensor. 
The resulting theory is called second order viscous hydrodynamics. One can shows that 
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for physically reasonable ranges of the second order coefficients the theory is causal 
In general there are large number of second order terms. A complete classification of the 
second order terms in a relativistic conformal fluid was recently given in |75j. Conformal 
symmetry implies that £ = and 5^T^ = —rja^. The second order correction is 



3 



(85) 

+ \ x a\a v ^ x + X 2 a { \Q^ x + X 3 Q { \Q u)x , 
where a^ u is the first order shear tensor defined above, 

AW = ^ a A^ [A aP + Ap a - ^A^A a ^j (86) 

denotes the transverse traceless part of A af3 and 

W = I A^A^ {d a u p - d p u a ) (87) 

is the vorticity. Equ. ( 1851) defines four new second order transport coefficients, tjj and 
Ai^^. These coefficient can be determined using kinetic theory [76] or the AdS/CFT 
correspondence [751 177] . 

Equation fl85l) is a constitutive relation that determines the stress tensor in terms 
of thermodynamic variables. Formally, we may replace time derivatives by spatial 
derivatives using the lower order equations of motion. Another option, inspired by 
the approach of Israel and Stewart [78], is to promote = 5T fll/ to a hydrodynamic 
variable. The equation of motion for is 



7T^ = - W ^ - Tjj 



+ — 7T A 7T ' 

This equation describes the relaxation of i\^ iv to the Navier-Stokes form —rja^. There 
are also a number of more phenomenological approaches that include some subset of 
higher order terms, for example the already mentioned Israel- Stewart formalism [7S] 
or the equations of Lindblom and Geroch [79], see [31] for a review. We note that 
whatever formalism is used, a necessary condition for the applicability of second order 
hydrodynamics is that higher order corrections are small, 6®T^ < 8^T^ < . 

Remarks: The second order formalism was initially developed for non-relativistic fluids 
by Burnett [80, 81J, see [82] for a review. Higher order hydrodynamic equations can be 
derived from kinetic theory by computing moments of the Boltzmann equation. This 
procedure is known as Grad's moment method |83j. It is not easy to find systems in 
which the second order theory provides a quantitative improvement over the Navier 
Stokes equation. An example is the work of Uhlenbeck, Foch and Ford on sound 
propagation in gases [EH [85]. Finally, we note that relativistic superfluid hydrodynamics 
was formulated by Carter, Khalatnikov and Lebedev [861 EZ] , see [88j [891 EH] for more 
recent studies that emphasizes the connection to effective field theory. 
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3.2. Diffusion 

An important diagnostic of the properties of a fluid is the diffusion of a dilute gas of 
impurities suspended in the fluid. We will see, in particular, that if the fluid is composed 
of quasi-particles then the diffusion of impurities and the shear viscosity, which is related 
to momentum diffusion, are closely linked. The two transport coefficients have the same 
dependence on the coupling constant, and their temperature dependence is the same up 
to kinematic factors. In a perfect fluid, however, this link may be broken: The diffusion 
constant goes to zero while the shear viscosity remains finite. 

We will assume that the number of impurity particles is conserved. The number 
density satisfies the continuity equation 

^ + V-j = 0. (89) 

If the number density varies smoothly then the current j can be expressed in terms 
of the thermodynamic variables. At leading order in derivatives of the density we can 
write j — — DVn, where D is the diffusion constant. Inserting this expression into the 
continuity equation gives the diffusion equation 

^ = DV"n. (90) 

A more microscopic view of diffusion is provided by studying the Brownian motion of 
an individual suspended particle. The motion is described by a stochastic (Langevin) 
equation 

^ = -vdp + m, mMf)) = - o- (9i) 

Here, p is the momentum of the particle, rjo is the drag coefficient, and £(£) is a stochastic 
force. The coefficient k is related to the mean square momentum change per unit time, 
3k, = ((Ap) 2 )/(At). The Langevin equation can be integrated to determine the mean 
squared momentum. In the long time limit (t ^> r/f) 1 ) the particle thermalizes and we 
expect that (p 2 ) = 3mT. This requirement leads to the Einstein relation 

Vd = tA^- (92) 

The relation between rjo and the diffusion constant can be determined from the mean 
square displacement. At late times ([Ax(t)] 2 ) = 6D\t\ and 

T 2T 2 

D = = . 93 

mr]D k 

A special case is the diffusion of large spherical particles suspended in a simple fluid. In 
this case the drag coefficient can be computed using the Navier-Stokes equation and the 
drag is related to the shear viscosity of the fluid, t]d = Gn^a/m, where a is the radius 
of the particles. This leads to a relation between the diffusion constant and the shear 
viscosity, D = T/{Qnrja). 
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3.3. Dynamic universality 

In the vicinity of a second order phase transition fluctuations of the order parameter 
relax slowly. This implies that order parameter fluctuations have to be included in 
the hydrodynamic description. The resulting hydrodynamic models describe universal 
features of transport phenomena near a continuous phase transition [91J. Dynamic 
universality classes, like the well-known static ones, depend on the symmetries of the 
order parameter and the number of dimensions. Universal aspects of transport also 
depend on the nature of the order parameter, whether it is conserved or not, and on the 
presence of couplings (non-vanishing Poisson brackets) between the order parameter and 
the conserved fields. In this section we will briefly review the hydrodynamic description 
of a simple fluid near the liquid-gas endpoint. This theory is known as model H in the 
classification of Hohenberg and Halperin |91j . We will see that critical fluctuations lead 
to a divergent shear and bulk viscosity at the liquid-gas endpoint. The hydrodynamic 
description of the superfluid-normal transition in liquid helium and dilute atomic gases 
is called model F. This model describes the divergence of the heat conductivity at the 
superfluid transition. 

Near the critical point sound modes (u ~ k) are higher in energy than diffusive 
modes (u ~ k 2 ), and the longitudinal components of the momentum density g can be 
neglected. A minimal model that describes the coupling of the order parameter <p to 
the transverse momentum density g^ is [9T] 



where P? = (5ij + VjVj/V 2 ) is a transverse projector, and ( gj are random forces, 
and the free energy JF T = T — is given by 



where h and A are external fields. The coefficients Ao and t]q are the bare values of 
the thermal conductivity and shear viscosity. Fluctuations cause the physical value of 
the zero frequency transport coefficients to diverge near the critical point. In order to 
study the critical behavior of the bulk viscosity the longitudinal component of g has to 
be included [52]. 

In a normal fluid the only conserved charges are the particle density, the energy 
density, and the momentum density. The order parameter is a suitable linear 
combination of the energy density and the particle density. In QCD the hydrodynamic 
variables include the chiral condensate, and the conserved energy density, baryon density, 
and isospin density. The QCD phase diagram is expected to have two critical points, 
one that corresponds to the endpoint of the nuclear liquid-gas phase transition, and 




(94) 



(95) 




(96) 



(97) 
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another one that is related to the endpoint of the first order chiral phase transition 
[93] . QCD hydrodynamics in the vicinity of the chiral critical point was analyzed by 
Son and Stephanov [91] who argue that the chiral endpoint, like the nuclear liquid-gas 
endpoint, is correctly described by model H. The values of the critical exponents can be 
determined using the epsilon expansion. The shear and bulk viscosity diverge with the 
correlation length £ as [91], E2] 

77 ~f*> 0,7-0.06), C~T C (x v ~2.8). (98) 

The critical endpoint is in the same static universality class as the Ising model and 
£ ~ r - 63 , where t — (T — T c )/T. We note that the diver gence in the bulk viscosity 
is much stronger than the divergence in the shear viscosity. These results demonstrate 
that, while there is empirical evidence for the suggestion that the viscosity minimum is 
located at endpoint of the liquid-gas phase transition (see Tabled]), this idea cannot be 
rigorously correct. Indeed, both rj/s and (/s diverge near the critical endpoint. 



3.4- Kubo relations and spectral functions 

Hydrodynamics is an effective description of the low energy, long wavelength response 
of a fluid. The transport coefficients appear as unknown constants in the hydrodynamic 
equations. These constants can be extracted from experiment, or computed from a more 
microscopic theory. The relationship between transport coefficients and correlation 
functions in a microscopic quantum field theory is provided by Kubo relations. We 
have seen that hydrodynamics fixes the low energy and low momentum behavior of the 
correlation functions of conserved charges, see equ. ( |62ll65l) . In the field theory these 
correlation functions can be computed using linear response theory. The response is 
governed by the retarded correlation function. In the case of shear viscosity the relevant 
correlation function is 

G% ' xy (oj, k) = -i fdtf d 3 a;e i( "'- k - x) 0(t)([T^(x,t),T^(O,O)]) , (99) 



where T^ v is the energy momentum tensor. The spectral function is defined by 
p(u>, k) = — 2ImCru(u;, k). The imaginary part of the retarded correlator is a measure 
of dissipation, while the correlation function S(cu, k) (see Sect. I3.1.ip is related to 
fluctuations. The relation between these two functions is called the fluctuation- 
dissipation theorem [95]. In the low frequency limit p(u;,k) = (u/T)S(uj, k). Matching 
the correlation function from linear response theory to the hydrodynamic correlator 
gives the Kubo relation 

V = lim lim £ K — -i . (100) 

w-*0 k^O 2uJ 

The formula for the bulk viscosity involves the trace of the energy momentum tensor 

C = - lim lim P - , (101) 

and analogous results can be derived for the thermal conductivity and diffusion 
constants. 
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T 


1.02 T c 1.24 T c 1.65 T c 


Tj/S 

C/s 


0.102(56) 0.134(33) 
0.73(3) 0.065(17) 0.008(7) 



Table 2. Lattice QCD results for the ratio of shear and bulk viscosity to entropy 
density in a pure gluon plasma. The calculations were performed for three different 
temperatures, given in units of the critical temperature T c . Data from |97[ I103j . 



The spectral function contains information about the physical excitations that 
carry the response. We will discuss this issue in more detail when we compare the 
strong coupling (AdS/CFT) and weak coupling spectral functions in Sec. HI Dispersion 
relations connect the spectral function to correlation functions with different analyticity 
properties. The Matsubara (imaginary energy) correlation function is 

GrfuJ (102) 

where u n = 27rnT is the Matsubara frequency. The imaginary time correlation function 
is given by 

Ge(t) = /^%rM,), (103) 
where the kernel K(uj,t) is defined by 

^ - ^™ = 11 + ^ ^ + ^ ■ <io4) 

and risiyj) is the Bose distribution function. Equation (I103j) is the basis of attempts 
to compute transport coefficients using imaginary time quantum Monte Carlo data 
[961 EH EH EH] • The idea is to compute Ge(t) numerically, invert the integral transform 
in equ. (11031) to obtain p(oo), and then extract transport coefficients from p'(0). The 
difficulty is that Ge{t) is typically only computed on a small number of points, and that 
the imaginary time correlator is not very sensitive to the slope of the spectral function 
at low energy. Many recent calculations make use of the maximum entropy method to 
obtain numerically stable spectral functions and reliable error estimates [1001 HOI] - It 
was also observed that one can minimize the contribution from continuum states to the 
imaginary time Green function by studying the correlators of conserved charges, energy 
and momentum density, at non-zero spatial momentum [1021 1103] . In more physical 
terms this means that one is extracting the viscosity from the sound pole rather than 
the shear pole. In Table [2] we summarize some recent lattice QCD results on the shear 
and bulk viscosity in the high temperature phase of pure gauge QCD. We observe that 
the shear viscosity to entropy density ratio is close to the conjectured bound l/(47r). 
The bulk viscosity is large in the vicinity of the phase transition but decreases quickly 
and becomes extremely small at T = 1.64T C . 
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3.5. Kinetic theory: Shear viscosity 

If the fluid can be described in terms of weakly interacting quasi-particles then the 
hydrodynamic variables can be written in terms of quasi-particle distribution functions 
/ p (x, t). In the case of a non-relativistic fluid the energy current, momentum current, 
and stress tensor are given by 

j;£(M) = J W EpVp4fp(?C,t) 1 (1 ° 5) 



, d^p 

9i(x,t) = I -—^mv Pji f p (x,t), (106) 



/d^p 
j^-^mv Pji v Ptj f p (x,t), (107) 

where Ep is the quasi-particle energy, and v Pt i = (dE p )/(dpi) is the quasi-particle 
velocity The equation of motion for / p (x, t) is the Boltzmann equation 

?A + v . Vfp + F . Vpfp = C [f p ], (108) 

where F is an external force and C[f p ] is the collision term. In local thermal equilibrium 
the distribution function is determined by the local temperature, chemical potential, 
and flow velocity We have 

#(x, t) = — — — - , (109) 

p exp((E p - v ■ p - n)/T) =F 1 

where the =F sign corresponds to bosons/fermions. Transport coefficients characterize 
how the distribution function relaxes to its equilibrium value if it is perturbed slightly 
away from it. We can write 

f p (x,t) = f°(x,t) + 5f p (x,t) (110) 

and linearize the Boltzmann equation in 5f p . In order to determine transport coefficients 
we also use a gradient expansion of the local velocity, temperature and chemical potential 
and linearize the Boltzmann equation in the "driving terms" VjU,-, VjT and V^/i. This 
procedure is known as the Chapman-Enskog method. In the next section we will describe 
the method in the case of phonon mediated transport in a superfluid, and then discuss 
some of the modifications that appear when studying high temperature Fermi and Bose 
gases as well as gauge theories. 

3.5.1. Phonons in dilute Fermi gases In the following we will concentrate on the shear 
viscosity of the low temperature, superfluid, phase of the dilute Fermi gas at unitarity. 
The calculation is similar to the computation of the shear viscosity of superfluid helium, 
but as explained in Sect. 12.31 the low energy effective theory of the dilute Fermi gas is 
more tightly constrained. We discuss the shear viscosity of liquid helium, as well as the 
viscosity of the low temperature (chiral symmetry broken) phase of QCD in Sect. 13.5.21 
The stress tensor of a phonon gas is 
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a) b) 




Figure 2. Leading order processes that contribute to the shear viscosity of the Fermi 
gas in the unitarity limit at low temperature (Fig. a) and high temperature (Fig. b). 
Dashed lines are phonon propagators and solid lines are fermion propagators. 



In order to study the shear viscosity we write Sf p = — + f p )/T with 

xip) = g(p) (piPj - \ 6 vP 2 ) ( v ^i + - l 6 n v ■ v ) • ( 112 ) 

Inserting this ansatz into the energy momentum tensor gives 

"=mf^^ 1+ ^- (113) 

The function g(p) is determined by the linearized Boltzmann equation. Linearizing the 
LHS of the Boltzmann equation in derivatives of v, fi, T gives 

dfp 2fp(^-~ >r fp) / 1 1 a \ 

—~ ~ cz— — - ' PijVa , (114) 

dt s 2E P T F3 J ' y 1 

where we have defined 

1 2 

Pij = PiPj - ^ijP 2 , Vij = ViVj + VjVi- -5 i:j V - v . (115) 

The RHS of the Boltzmann equation contains the collision term C[f p }. In the present 

case the dominant contribution arises from binary 2 «-> 2 collisions. The linearized 
collision term is 

CW/ P ] * / dT{k-k\p'){l + f k )f kl f pl \M\ 2 (116) 



x 



g(p)Pij + gitykij - g{k')k' i:j - g(p')p'i 



= 2ET P Cij\g(p)]vij, 

where M. is the scattering matrix element, 

dT(k;k\p')=( I] J^S^r) ^W5^\p + k~k'-p') (117) 
\q=k,k>, P > \ Z71 ) J 

is the phase space, and we have defined the linearized collision operator Cij[g(p)]. The 
linearized Boltzmann equation can now be written as 

C-Mp)] = yP^ (118) 
This result can be used to rewrite the relation for the viscosity in equ. (11131) as 
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The two relations equ. (1113)) and ()119p can be used to derive a variational estimate 
of the shear viscosity. We can view equ. (1113)) as an inner product with measure 
+ f°)/(2E p ) and write 

V=\{X\g) ) (120) 

where X = (c 2 /T)pij and g = g(p)pij. The linearized Boltzmann equation is C\g) = \X) 
and equ. ( 11191) can be written as r/ = | (g\C\g). The linearized collision operator C is a 
hermitean, negative semi-definite operator. The zero eigenvalues of C correspond to the 
conservation laws for energy, momentum, and particle number. Consider a variational 
ansatz \g var ) for the exact solution \g) of the linearized Boltzmann equation. The triangle 
equality implies 

(9var\C\g var )(g\C\g) > (g va r\C\g) 2 = (g var \X) 2 . (121) 
Using r] = | (g\C\g) we get 

2 (g va r\X) 2 , , 

V -5 (a cJ (122) 

° \gvar\^ \gvarj 

This result is, of course, not a lower bound on the exact value of rj, but it is a bound 
within the approximation that is used to compute the collision term. A popular choice 
for g var is the driving term X. This ansatz provides a good estimate in the case of 
non-relativistic particles interacting via short range interactions, as well as for gauge 
boson exchanges in QCD, but not in the case of phonon scattering^]. 

A systematic method for improving the variational estimate is based on orthogonal 
polynomials. We can construct a complete set of polynomials that are orthogonal with 
respect to the inner product defined in equ. (11201) . In non-relativistic physics these 
polynomials are known as Sonine polynomials [10TJ and suitable generalizations can 
been constructed for Bose and Fermi gases |108] . We now fix an integer N and expand 
the solution of the linearized Boltzmann equation in the first N polynomials. At finite N 
solving the Boltzmann equation reduces to the problem of inverting an N x N matrix. 
The solution provides a variational estimate for rj which becomes exact as N — > oo. 
Convergence is usually quite fast. 

To complete the calculation of the shear viscosity we need to compute the scattering 
amplitude A4. The collision operator at leading order in T//x is determined by the 
scattering amplitude at leading order in where q = p,p',k,k'. The amplitude 
is given by the diagrams in Fig. [2|i with vertices and propagators determined by the 
effective lagrangian given in equ. (1281) . The expression for M. is not very instructive and 
can be found in [108j . The best estimate for rj is obtained by using g(p) ~ p~ x . We find 

rj = 9.3 x 10- 6 ^ , (123) 
cj T 5 

II A detailed discussion of upper and lower bounds on transport coefficients can be found in [10411105) . 
We also refer the reader to comparisons of the variational results with exact solutions of the Boltzmann 
equation [106] . 
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where £ ~ 0.4 is the universal parameter we introduced in equ. 0231) . In the low 
temperature limit the entropy density is dominated by the phonon contribution 
2tt 2 T 3 

S = -7T-T- ( 124 ) 
45 c 3 s 

The ratio rj/s drops sharply with temperature. Extrapolating to T = T c ~ 0.15Tp gives 
rj/s ~ 0.8, with very large uncertainties. 

3.5.2. Phonons and rotons in liquid helium, pions in QCD The calculations of shear 
viscosity of liquid 4 He below the A point is similar to the computation of rj in the 
superfluid Fermi gas. The main difference is that close to T c it is important to include 
the roton contribution. Rotons form a dilute gas, and unlike phonons, their cross section 
is approximately constant. As a consequence the roton viscosity is independent of 
the roton density, see the discussion below equ. (j2j). The typical roton momentum 
is determined by the roton minimum of the dispersion relation and depends only 
weakly on temperature. This implies that the roton viscosity is almost temperature 
independent. The value of the roton viscosity depend on the poorly known roton-roton 
interaction. A fit to experimental data for the shear viscosity below the lambda point 
gives rj r ~ 1.2 ■ 10~ 5 poise. The leading correction to the roton term comes from phonon- 
roton scattering. Landau and Khalatnikov find [25] 

where A is the roton energy defined in equ. (|I6]) . A is a constant, and Q/Q p h is the ratio 
of the roton-roton and roton-phonon relaxation rates. This ratio scales as T 4,5 exp(A/T). 
For T < 0.9 K we can use 6 ^> Q p h and the temperature dependence of the phonon- 
roton term is governed by the T~ 05 exp(A/T) term. For T < 0.7 K phonon-phonon 
scattering dominates and the viscosity scales as T~ 5 , as in the previous section. At 
even smaller temperature, T < 0.5 K, phonon splitting, also known as Beliaev damping, 
becomes important and the temperature dependence changes to rj ~ T _1 |109j . The 
roton contribution to the entropy density is 

»= (W)W% ( 1 + 2AM-r)' < 126 > 

where m* and po are given in equ. (Tl6l) . The phonon contribution is given by equ. (11241) 
with c s = 238 m/sec. If we push equ. fl 125[) and (11261) to the limit of their applicability, 
T ~ 2 K, we find 77/s ~ 2. 

The computation of the shear viscosity in low temperature QCD also proceeds 
along similar lines. The analog of the phonon in QCD is the pion, and pion interactions 
are governed by the effective lagrangian given in equ. fT46l) . The pion is not massless, 
m w = 139 MeV. At very low temperature, T m^, the pion scattering amplitude 
is approximately constant and the viscosity is only weakly temperature dependent. At 
higher temperature we can set m n ~ and the scattering amplitude is energy dependent. 
The main difference as compared to phonon scattering is that the four-pion interaction 
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is of the form (<j)d(j)) 2 instead of (<90) 4 , and that there is no three-pion interaction. As a 
consequence the tot scattering matrix element scales as the second power of the external 
momenta. The pion entropy is given by equ. ( 11241) with c s = cj \/3 and an isospin 
degeneracy factor 3. An approximate calculation of the ratio r\js gives [1101 1111] 

v = (127) 

s 16ttT 4 v ; 

Variational solutions of the Boltzmann equation reported in [1 12] give rj/s ratios that 
are about five times larger. The first study of the shear viscosity of a pion gas can be 
found in [113j . More detailed investigations of the viscosity of hadronic mixtures were 
published in [rnJlfTTi] . 

3.5.3. Non-relativistic atoms: Dilute Fermi gases and 4 He The shear viscosity of the 
dilute Fermi gas at high temperature is determined by binary scattering between the 
atoms. The s-wave scattering matrix is 

47T 1 

M = — —. -, (128) 

m 1/a + iq 

where q is the relative momentum. In the unitarity limit a — > oo the scattering amplitude 
diverges as 1/q in the limit of small momenta. For T ^> Tp the infrared divergence 
is effectively cut off by the thermal momentum (mT) 1 / 2 . The viscosity in the high 
temperature limit is |115j 1 1 6J 

1~Kf. (129) 

This result is based on the variational function g(p) ~ 1. Corrections due to more 
complicated distribution functions are small, A77/77 < 2% |117j . The scaling of rj with 
temperature can be understood as a combination of the T 1//2 scaling of a dilute hard 
sphere gas (see Sect. [1]) with an extra factor (mT) from the 1/g 2 behavior of |.M| 2 . The 
high temperature limit of the entropy density is that of a classical gas 



2V2 
3^ 



{mT F fl 2 



\ /3v/?T 3 / 2 \ 5 



(130) 



Combining equ. (11291) and (I130p gives r]/s ~ x 3 / 2 /log(x) with x = T/T F . The classical 
expression for the entropy becomes unphysical (negative) for T ~ T c . Extrapolating to 
T ~ 2T C gives t]/ s ~ 0.5. 

The shear viscosity of helium is governed by scattering in the potential given in 
equ. (fl~5l) . In the high temperature limit the dominant contribution does not come from 
the Van der Waals tail, but from the repulsive short range contribution. For a potential 
of the form V ~ r~ v the viscosity scales as T s with s = | + [118] . In the case of 
a Lennard- Jones (6-12) potential this implies i] ~ y - 68 . A somewhat better fit to the 
data is provided by 



T 



0.647 



V = Vo[7 F ) (131) 
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Figure 3. Leading order processes that contribute to the shear viscosity of a pure gluon 
plasma. The coefficient k defined in equ. (|133[) is determined by the i-channel diagram. 
The full leading order result, including the coefficient fj,*, requires the remaining 
diagrams, as well as gluon bremsstrahlung from the external legs (not shown). 



with 7/0 = 1.88 ■ 10~ 5 Pa-s and T = 273.15 K. The entropy density is given by the 
classical result, equ. (11301) . For T = 10 K we get rj/s ~ 4, and extrapolating to T = 4 
K gives rj/s ~ 1.5. Very accurate calculations that are based on realistic potentials and 
include higher order terms in the density can be found in j!19j . These calculations are 
reliable down to about 10 K. 



3.5.4- Gauge fields in QCD The shear viscosity of a quark gluon plasma is determined 
by binary quark and gluon scattering. We first consider a pure gluon plasma. The 
leading order gluon-gluon scattering diagrams are shown in Fig. HI The squared tree 
level amplitude is 

9g 4 / ut us ts' 

where g is the gauge coupling and s, t, u are the Mandelstam variables. The differential 
cross section diverges for small momentum transfer q as 1/g 4 . This is the standard 
Rutherford behavior, which arises from t-channel gluon exchange. In the calculation 
of the shear viscosity the differential cross section is weighted by an extra factor of 
(1 — cos 9), where 9 is the scattering angle. The quantity o~t = J dcos9 (do~)/(d cos 9){\ — 
cos 9) is known as the transport cross section. The transport cross section diverges 
logarithmically at small 9. This divergence is regulated by medium corrections to the 
gluon propagator, see equ. (I39ll40p . Electric gluon exchanges are screened at a distance 
To ~ m^ 1 , and the electric contribution to o~t is proportional to (/ 4 log(m/j). There is 
no static magnetic screening, but gluons with energy uj are dynamically screened at a 
distance r ~ (cum^,) -1 / 3 . After integrating over energy the magnetic contribution also 
scales as g 4 log(m/j). Combining electric and magnetic t-channel exchanges gives [16|[T7] 

T 3 

V = k— — rT7 r, (133) 

g 4 log{p,*/m D ) 

where k = 27.13. We will specify the coefficient fi* below. This result corresponds to an 
optimized trial function xip) = ^ip)Vij v iji but the simple approximation A(p) ~ const 
agrees with the exact result to better than 1%. In order to compute the shear viscosity 
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of a quark gluon plasma we have to include t-channel diagrams for quark-quark and 
quark-gluon scattering. The result is of the same form as equ. H133[) with [T7j 

k(N f ) = (27.13,60.81,86.47,106.67), (^ = 0,1,2,3). (134) 

Note that k increases with Nf faster than the total number of degrees of freedom. This is 
related to the fact that quarks have smaller color charges than gluons, which implies that 
quark-gluon scattering amplitudes are suppressed relative to gluon-gluon amplitudes. 

In order to make an absolute prediction for the shear viscosity we need to determine 
the constant fi* in equ. (11331) . This coefficient receives contributions from s and u- 
channel gluon exchanges. These contributions are straightforward to include. A more 
difficult problem arises from the fact that fi* is sensitive to soft (q ~ mo) binary 2^2 
scattering followed by collinear 1 — > 2 splitting. The inverse mean free time for this 
process is given by r _1 ~ g 4 T 3 /m 2 D x g 2 ~ g 4 T, comparable to the transport mean 
free time r^. 1 ~ T 4 /r] ~ g 4 T. Since the scattering angle is zero collinear splitting does 
not directly contribute to shear viscosity, but it degrades the momentum and assists in 
randomizing the momentum distribution in subsequent binary collisions. 

The difficulty with collinear splitting is that the formation time of the emitted gluon 
is of order l/(g 2 T). This is the same order of magnitude as the quasi-particle life time 
given in equ. (14*31) . which implies that kinetic theory is breaking down. Arnold, Moore, 
and Yaffe showed that if interference between subsequent gluon emission processes, the 
Landau- Pomeranchuk effect, is taken into account an effective Boltzmann equation with 
2^2 and 1 — > 2 collision terms can be derived pL20j. Arnold at al. find [18] 

fj,*(N f = 0) = 2.765 T. (135) 

They also show that fi* is only weakly dependent on the number of flavors, fi*(Nf = 
3) = 2.957 T, and compute additional terms in an expansion in inverse logarithms of 
fi*/m D . 

The entropy density of the quark gluon plasma is given by 

s=^- (2(N 2 -l) + 7 -W f ^ T 3 . (136) 

Higher order corrections to the entropy density are large, but the situation in the regime 
T > 2T C can be improved using resummation schemes, see Fig.[TJ The resummed entropy 
differs from the free gas result by no more than 15% for T > 2T C . The magnitude of 
higher order corrections to the viscosity is not known, but next-to-leading order results 
for the heavy quark diffusion constant suggest that higher corrections to transport 
coefficients are large |121j . 

The leading order QCD result is shown in Fig. HI Clearly, t]/s is strongly dependent 
on the coupling, and without performing higher order calculations it is not clear what 
value of a s one should use at a given temperature. An interesting perspective is provided 
by exact results for rj/s in the strong coupling limit of Af = 4 SUSY Yang-Mills theory, 
see Sect. HI These results can be compared to weak coupling calculations based on 
kinetic theory |123j . The weak coupling result for rj/s in the M = 4 theory is smaller 
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Coupling a 3 Coupling A=N c g 2 

Figure 4. Shear and bulk viscosity to entropy density ratio in QCD (left panel) and 
Af = 4 supersymmctric Yang-Mills theory (right panel) . The left panel shows the shear 
and bulk viscosity to entropy density ratio in QCD with Nf = 3 flavors as a function 
of the strong coupling constant a s , from |122j . The right panel shows the ratio rj/s in 
TV = 4 SUSY Yang Mills theory as a function of the 't Hooft coupling A = g 2 N c . The 
solid line shows the weak coupling result, the dotted line is an extrapolation of the 
weak coupling result to the strong coupling regime, the dashed lined is A — * oo result 
from the AdS/CFT correspondence, and dash-dotted line is the leading correction to 
the strong coupling result, from [123] . 

than the corresponding ratio in QCD by a factor ~ 1/7. This is related to the fact that 
in the M = 4 theory all states are in the adjoint representation, and that the theory 
contains extra scalars. Both of these differences lead to larger cross sections. 

Weak and strong-coupling results for rj/s as a function of the 't Hooft coupling 
A = g 2 N c in SUSY Yang-Mills theory are shown in the right panel of Fig. HJ We observe 
that r]/s in the Af = 4 theory reaches the strong coupling limit when extrapolated to 
a 't Hooft coupling A = g 2 N c ~ 12. As discussed in Sec. 12.51 this is a large N c result. 
Naively extrapolating to N c = 3 the value A ~ 12 corresponds to a s = # 2 /(47r) ~ 0.3. 
We also note that the value of 't Hooft coupling at which the weak coupling result for 
r]/s reaches the strong coupling limit is larger than the coupling A ~ 5 at which the 
corresponding expression for the entropy reaches the strong coupling limit s/sq = 0.75, 
see Fig. [U If we consider s/sq = 0.8 to be the "QCD-like" point, then we should restrict 
ourselves to A < 5. In this case rj/s does not drop below 0.5. 

3.6. Kinetic theory: Other transport properties 

3.6.1. Bulk viscosity Bulk viscosity measures the amount of energy dissipated as a 
fluid is slowly expanded or compressed. In a conformally invariant system changing all 
the momenta and positions by a constant scale factor connects equilibrium states and 
the bulk viscosity must vanish. In kinetic theory bulk viscosity is typically sensitive to 
processes that change the particle number or the composition of the system. The kinetic 
theory prediction for bulk viscosity is proportional to the corresponding relaxation time, 
and to deviations from conformality in the equation of state. Depending on the interplay 
between these two effects, the temperature dependence of the bulk viscosity can differ 
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dramatically between different fluids, and between shear and bulk viscosity. There 
are many fluids for which bulk viscosity is not an important source of dissipation, 
either because they are approximately incompressible, like water, or because the fluid is 
compressible but approximately scale invariant, like the QGP plasma. On the other, we 
have seen that bulk viscosity is the dominant source of dissipation near a second order 
phase transition, see equ. fl98|) . 

The Fermi gas at unitarity is exactly conformal and the bulk viscosity in the 
normal phase vanishes. In the low temperature phase conformal invariance requires 
Ci — C2 = 0, but (3 can be non-zero [124] . This coefficient was recently computed in 
|125j . The result is sensitive to non-linearities in the phonon dispersion relation. If 
1 — > 2 phonon splitting is kinematically allowed then (3 ~ T 3 , where the constant of 
proportionality depends on the curvature of the dispersion relation. The bulk viscosity 
of liquid helium was calculated by Khalatnikov [25]. As in the case of shear viscosity 
the main contribution comes from phonons and rotons. Khalatnikov finds that (2, 
the bulk viscosity of the normal component, is about an order of magnitude bigger 
than rj. The other two bulk viscosities, £1 and £3, involve motion of the normal fluid 
relative to the superfluid. They have different physical units, and cannot be directly 
compared to £2- The linear combination that appears in the damping of second sound 
is = (2 + R 2 Cz ~~ 2pCi- At normal density there are significant cancellations between 
these terms and (2 ~ (p 2 C3 — 2pCi)- The bulk viscosity of helium vapor is small. Note 
that the bulk viscosity of most gases is dominated by internal excitations, like rotational 
and vibrational modes. 

The bulk viscosity of a pion gas at low temperature was computed by Chen and 
Wang [126]. They find that the bulk viscosity scales as ( ~ T 7 / f£ (up to logarithms). 
The bulk viscosity of the high temperature QGP phase was calculated by Arnold, Dogan, 
and Moore [122] . The result is 

C = ; 7 "f , , (13T) 

where A = 0.443 and p* = 7.14 T in pure gauge QCD. In full QCD with N f = 3 quark 
flavors A = 0.657 and p* = 7.77 T. We observe that ( scales as a 4 x rj. The trace 
anomaly e — 3P is proportional to ot 2 s , so bulk viscosity scales like the shear viscosity 
times the second power of the deviation from conformality. This is in agreement with 
an a simple formula proposed by Weinberg [127] . ( ~ (cf — |) 2 ?7- However, Weinberg's 
relation is known to be violated in some theories, see |128j for an example. 

3.6.2. Diffusion The diffusion of of impurities in liquid helium has been studied in 
some detail. Of particular interest is the behavior of dilute solutions of 3 He in 4 He. 
At low temperature the diffusion constant is determined by scattering off phonons and 
D ~ 1/T 7 [25], 1129] . At high temperature diffusion is governed by scattering between 
atoms and D ~ T 1+s with s = \ + for a 1/r" potential |118j . We conclude that 
the temperature dependence of the diffusion constant is identical to that of the shear 
viscosity. In case of the unitary Fermi gas one can make use of the fact that the number 
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of spin up and spin down fermions is separately conserved, and study the diffusion of 
minority spin down particles in a background of majority spin up fermions [130j . This 
process contains important information about the interaction between the different spin 
states, but it is not directly related to the viscosity of the spin balanced gas. 

The diffusion constant for heavy quarks in a quark gluon plasma can be determined 
by computing the mean square momentum transfer per unit time, see Sec. 13. 21 For 
approximately thermal heavy quarks the diffusion constant is dominated by heavy quark 
scattering on light quarks and gluons, qQ — > qQ and gQ — > gQ. As in the case of shear 
viscosity the most important Feynman diagrams involve t-channel gluon exchanges. 
Since the heavy quark is slow the dominant interaction is electric gluon exchange and 
the cross section is regularized by Debye screening. The leading order result is [131[ 1132] 



D 367F 
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C F g*T 

where C F = (N 2 C - 1)/(2JV C ) and c = 0.5 - ^ E + C'( 2 )/C(2). We note that the 
diffusion constant has the same parametric dependence on the coupling as the shear 
viscosity. The relaxation time rj^ 1 scales as Mg/(T 2 g 4 log(g), which is larger by a factor 
Mq/T compared to the hydrodynamic relaxation time r]/(sT) ~ 1 / (T g A log(g)) . This is 
confirmed by numerical estimates, which give — 6.7/T ~ 7 fm for charm quarks at 
T = 200 MeV [132] . For comparison, the hydrodynamic relaxation time is r)/(sT) ~ 1 
fm (for r]/s ~ 1). 

3.6.3. Thermal conductivity Thermal transport in superffuid helium is a complicated 
process. In a superfluid heat transport can take place by a process similar to internal 
convection where the superfluid moves relative to the normal fluid. Only the normal fluid 
carries entropy and as a result heat is carried along with the normal component. The 
convective contribution to heat flow is controlled by the shear viscosity of the normal 
fluid. Within the normal fluid heat is carried by phonons and rotons. Khalatnikov 
showed that there is no heat transport in a gas of phonons with exactly linear dispersion 
relation [25]. The thermal conductivity of the normal fluid is dominated by rotons and 
phonon-roton scattering. This situation is somewhat similar to heat transport in a solid. 
At very low temperature, heat transport is ballistic and the entropy is carried by a net 
flow of phonons along the temperature gradient. At higher temperature non-linearities 
in the phonon dispersion relation and the "umklapp" process play a role. 

The situation at high temperature is much simpler. Heat flow is a diffusive process, 
and the thermal conductivity is determined by scattering between atoms. A simple 
mean free path estimate analogous to equ. (j2j) is 

K = \nc p pl mfp , (139) 

where c p is the specific heat at constant pressure. This estimate suggests that the ratio 
of the shear and thermal diffusion constants, the Prandtl number 

Pr = ^, (140) 

K 
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is close to one. At large T the thermal conductivity of helium scales as T s with 
s = | + ^y, as in the case of shear viscosity. The Prandtl number is approximately 
constant, Pr ~ 2.5. 

Most studies of the thermal conductivity of a quark gluon plasma have focused 
on the regime of very high baryon density. In the limit jj, ^> T, where \i is the quark 
chemical potential, the thermal conductivity scales as k ~ /i 2 /a 2 . |133j . In the opposite 
limit T ^> jj, there is an old relaxation time estimate k ~ T 4 / (a 2 /i 2 ) [9j. Note that while 
k diverges as \i — > 0, the dissipative contribution to the baryon current, equ. ( JHTj) . is 
finite. 

4. Holography 

In kinetic theory conserved charges are carried by well defined quasi-particles. The 
time between collisions is long compared to the quantum mechanical scale, h/T, and 
quantum mechanical interference between scattering events is not important. In the 
strong coupling limit quantum mechanical effects are large and quasi-particles lose their 
identity. A powerful new tool to study transport phenomena in this regime is the 
AdS/CFT correspondence [201 flMj 1155] . 

The AdS/CFT correspondence is referred to as a holographic duality - it relates 
string theory on a certain higher dimensional manifold to four dimensional gauge theory 
on the boundary of this space. The correspondence is simplest if the field theory is 
strongly coupled. In this limit the string theory reduces to a classical gravitational 
theory. The holographic correspondence then implies that a four dimensional field theory 
is capable of encoding gravity in five dimensions. The idea of a correspondence between 
field theories and higher dimensional gravity originated from developments within string 
theory, but there are precursors to the correspondence that come from the physics of 
black holes. It has been known for some time that black holes carry entropy, and that 
the entropy is proportional to the area, and not the volume of the black hole. It was 
also known that the evolution of black holes respects the second law of thermodynamics, 
and that it can be described by treating the event horizon as a physical membrane with 
well defined transport properties like electric conductivity and shear viscosity [136J. 

The best studied example of the AdS/CFT correspondence is the equivalence 
between M = 4 Super Yang Mills theory (see Sect. 12.51) and string theory on AdS$ x S5. 
For our purposes the dynamics only involves AdSs. This is a 5-dimensional space, which 
in AdS / CFT terminology is called the bulk. The dual field theory exists on the boundary 
of this space, which is 3+1 dimensional Minkowski space. The gauge gravity duality 
works as follows: Classical gravity equations of motion are solved in the 4+1 dimensional 
curved geometry of AdSs. Fluctuations of gravitational fields in the bulk induces charges 
on the 3+1 dimensional boundary. The dynamics of 3+1 dimensional boundary theory is 
the strongly coupled conformal field theory which we wish to study. Transport properties 
of the boundary theory can be determined by perturbing the boundary charges with an 
external field which then propagates into the bulk. The response of the induced charges 
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to the applied field determines the transport coefficients. For each conserved charge 
of the field theory there is a corresponding field in the gravitational theory. The field 
corresponding to the stress tensor T^ v is the graviton h) iV ', and the field corresponding 
to the conserved R charge current is the five dimensional Maxwell field A 11 . 

The AdS/CFT setup is analogous to a parallel plate capacitor. Electromagnetic 
fields in the bulk, the space between the plates, induce surface charges on the boundary. 
Fluctuations of the bulk field create fluctuations of the surface charges, and correlation 
functions of the surface charges can be related to normal modes of the bulk field. What 
is remarkable about the AdS/CFT correspondence is that the gravitational theory 
in the bulk defines a local field theory on the boundary, and that there are classical 
gravitational field configurations that correspond to field theories at finite temperature. 
These configurations can be used to study dissipative phenomena in the boundary field 
theory. 

The gravitational field configuration relevant to field theories at finite temperature 
is an AdSs black hole. In the black hole geometry the gravitational field is non zero 
as we approach the boundary of the 4 + 1 dimensional space. This gravitational field 
is balanced by a non-zero stress tensor in the boundary field theory - the gravitational 
setup corresponds to the dynamics of a field theory with a non- zero density matrix. The 
event horizon of the back hole spans three spatial dimensions in the bulk and radiates 
at the Hawking temperature T#. The black hole fills AdSs with a bath of gravitational 
radiation, and the temperature of the heat bath is identified with the temperature of 
the boundary field theory. The dynamics of graviton propagation in the black hole 
background determines stress tensor correlators at finite temperature in the boundary 
field theory. These correlators determine the shear viscosity according to Kubo formulas. 

There is a vast amount of literature on the AdS/CFT correspondence. A detailed 
review with extensive references is |137j . and more pedagogical reviews can be found in 
[138[ 11391 1140] . Reviews with an emphasis on transport phenomena are [2E1 EH] . Here we 
will concentrate on a few selected issues that are relevant to this review. First, we will 
explain the calculation of the shear viscosity and the spectral weights of strongly coupled 
fluids. Then we will comment on the conjectured viscosity bound, and the calculation 
of other transport properties. Finally, we will review the derivation of higher order 
fluid dynamics using holography, and summarize some recent attempts to extend the 
correspondence to non-relativistic theories. 

4-1. The equation of state from holography 

AdS 5 x S 5 is the product of five dimensional Anti-deSitter Space (AdS) and a five- 
sphere. Anti-deSitter space is a simple solution of the source free Einstein equation 
with a negative cosmo logical constant. Note that, on large scales, our universe is a 
approximately a four dimensional deSitter space. The geometry of AdS 5 x S 5 is described 
by the metric 




(141) 
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Here, dQf is the metric of the five-sphere, and (t, x, r) are the coordinates on AdSs. The 
coordinate r is referred to as the "radial" AdSs coordinate. The limiting value r — > oo 
is the "boundary" of AdSs. A fixed r slice of AdSs is a 3 + 1 dimensional flat Minkowski 
space, but the five dimensional space is curved, with a constant negative curvature. 
L is the corresponding curvature radius. We require that L is large compared to the 
string length £ s which guarantees the validity of the classical approximation. In the 
AdS/CFT correspondence L is related to the coupling constant of the M = 4 gauge 
theory, A = gy M N c , through the relation (L/£ s ) 4 = A. The classical approximation to 
the gravitational theory is reliable if the field theory is strongly coupled. The classical 
fields can be expanded in S5 spherical harmonics. At strong coupling higher harmonics 
are separated by a large gap, and we will ignore the S 5 from now on. 
The metric of an AdSs black hole is 

,2 . L 2 

L 2 v J v — y 1 f^y 2 

where /(r) = 1 — (r /r) 4 . The black hole horizon is a 3 + 1 dimensional surface at 
r = vq. The horizon radius is related to the Hawking temperature of the black hole by 
TqTi/L 2 = TrTff. This formula is an example of a general radius-energy relation in the 
AdS/CFT correspondence. A modification of the AdS geometry at radius r corresponds 
to a modification of the field theory at an energy scale rh/L 2 . It is convenient to perform 
a change variables u = (r /r) 2 and write the metric as 

^TLf ( ,, AM3 , ^ L 2 



ds " = T2 [-f( r )dt 2 + rfx 2 ) + T7 -^- 2 dr 2 , (142) 



ds 2 



u (-mdt' + d^ + ^j^du 2 , (143) 

where f(u) = 1 — u 2 . Now the horizon is at u — 1. The boundary limit is found by 
evaluating all quantities at u — e and then taking the boundary limit e — > 0. 

As discussed in the introduction to this section, the modified metric implies that 
there is an induced stress tensor at the boundary, u = e. This is an important point, 
and we will compute the induced stress tensor in two different ways. First, we will 
determine it by varying the action with respect to the boundary metric. This is the 
standard method by means of which one can determine the source of a given gravitational 
field. The only unusual ingredient is the fact that induced stress tensor is located on 
the boundary. We will provide an alternative derivation based on the analogy with the 
induced surface charge in electrodynamics below. 

The boundary metric is related to the metric of the five dimensional theory G^ u 
by the AdS scale factor 

9 ^ = J^fl^ G ^- (M4) 

Here and below Greek letters denote four dimensional indices (x^ 1 ) = (t, x, y, z) while 
Roman letters denote five dimensional indices (x M ) = (a; M ,-u). Near the boundary the 
metric can be written 

g^ = g° u + u 2 B^ + 0(u 4 ), (145) 
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where g° is interpreted as the metric of Af = 4 gauge theory. Usually g° is simply 
r/p. We will see that the coefficient of u 2 determines the induced stress tensor on the 
boundary. 

The induced stress tensor is 

-2 5S 



{Tp,,) = lim 



(146) 



+0 V-gSg^ 

where \f—g = (—detg^) 1 ^ 2 . The action is a sum of the Einstein-Hilbert action, the 
Gibbons-Hawking- York boundary term, and counter terms which are needed to render 
the action finite in the limit u — > 0, 

S = S EH + S GH + S CT . (147) 

The Einstein-Hilbert action is 

Seh = ^[ d 5 x^(n + 2A) , (148) 

2^5 Jm 

where 1Z is the Ricci scalar and A = 6/L 2 is the cosmological constant. The five 
dimensional Newton constant is related to the number of colors in the field theory, 
= N 2 /(4tt 2 L 3 ). The Gibbons-Hawking- York j!41l I142j boundary action is 



S gh = 7t^ d i x^2K, (149) 

2k§ JdM 

where we have defined the boundary metric 

7^ = <^U> (150) 

and K is the trace of the extrinsic curvatur^l The boundary term guarantees that 
the variation of the action with respect to the 5-dimensional metric gives the Einstein 
equations in the bulk provided the variation vanishes on the boundary. Without the 
boundary action one also has to require that derivatives of the variation vanish on the 
boundary, see [143J. Finally, the counter term 

ScT = -jf d 4 x^, (151) 

L JdM 

is needed to render the action finite in the limit u — > 0. Notice that the counter term 
is independent of temperature. With these definitions, the variation relates the stress 
tensor to the extrinsic curvature 



(7^) = — -J Hm 



1 i^L) 



2 



3 



(152) 



Kg u->0 u 

Substituting the black hole metric equ. (11431) and using the definition of the extrinsic 
curvature we have 



e N 2 

(Tp,) = diag(e, p,p,p) , -= p = — |(ttT). (153) 

We find that e = 3P, as expected for a scale invariant theory. We can compare the 
coefficient of T 4 to its value in the non-interacting theory. J\f = 4 SUSY QCD has 

» More explicitly, K = G^V^n,/ with n M an outward directed normal to the boundary of the AdS 
space, n M = -VG^S 5M . Note that = V M ?v = -n„r^ = n u d u G^. 
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8 (AT? — 1) ~ 8Af? bosonic and fermionic degrees of freedom. The contribution of a 
massless fermion to the pressure is 7/8 of that of a massless boson, see equ. fl35|) . 
Equation (11531) shows that the pressure in strongly coupled N = 4 SUSY QCD is three 
quarters of the Stefan Boltzmann value. 

We can also obtain equ. (1152ft in analogy with the induced surface charge on a 
capacitor plate. Consider a plate that spans the x-y plane. The surface charge density 
is related to the jump of electric field across the plate 



a 



E z 



(154) 

where [E z ] = — E z _ notates the jump. The analogous formulas in the gravitational 
theory are known as junction conditions |144j . Integrating the Einstein equations across 
a Gaussian pill box relates the surface stress r M „ to the jump in extrinsic curvature 



(155) 



Thus the particular combination of extrinsic curvature plays an analogous role to the 
normal electric field, i.e. a combination of —K^ v = n u T^ u is the analog of n • E. If we 
have a semi-infinite metal block with surface charge density a, then the outgoing electric 
field is related to the surface charge E z = a. By analogy, we associate the outgoing flux 
of extrinsic curvature at u = e with the stress tensor in the gauge theory 

= -- 2 ^ (K» v - K5%) . (156) 

Then taking the boundary limit u — » 0, we tentatively define the stress 

n = ~- 2 Km (K\ - K5» v ) . (157) 

Substituting the black hole AdS metric into this expression gives a divergent result. 
Nevertheless, the difference between this stress and the stress determined with the 
vacuum AdS metric equ. (11411) is finite 

< r " > - <^L,»„, = ^ (*". - K". + I s ") ■ < 158 > 

After lowering the indices of K^ v with 7^ = [(nTL) 2 /u] 77^, this equation is the same 
as derived previously in equ. (I152p . 



4-2. Shear viscosity from holography 

In the previous section we computed the average stress tensor on the boundary, 
(T^x, £)), which is a one point function of the conformal field theory. By Kubo's 
formula, equ. (llOOp . the shear viscosity can be related to a retarded two-point function. 
We will determine this function using linear response theory. Momentarily ignore the 
fifth dimension and consider turning on a time varying gravitational field h° (u) in 
the usual four dimensional field theory. This time varying gravitational field induces a 
deviation from the equilibrium stress tensor in the same way that a time varying electric 
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field induces a net current. According to linear response theory, the expectation value 
of the stress energy tensor is 

{T xy {u)) Ky = T%{u) + G R (u) h° xy {u) , (159) 

where — (e + p)u x u y + pg xy = ph° xy (uj) is the equilibrium stress tensor, and Gr{uj) 
is the equilibrium retarded correlator defined in equ. (199]) . Kubo's formula dictates the 
functional form of this correlator in the small frequency limit, Gr{uj) = —ioorj. Thus the 
average stress tensor in the presence a time varying gravitational field is 

(T xy (u)) Ky = ph° xy - iur) h° xy (u) . (160) 

Now consider the small fluctuations of the metric field H xy (u,u) around the black 
hole metric equ. f!143j) of the five dimensional theory. The equation of motion for the 
gravitational fluctuation is found by linearizing the Einstein equations 

K MN -^G MN (K + 2A) = 0. (161) 

After a modest amount of algebra, the lZ xy equation becomes an equation for h xy = 
uH xy (u,u)/(irTL) 2 

1 + u 2 u 2 
K y (uj, u) - —^h' xy (u, u) + — — ^ /^(fa,, u ) = 0, (162) 

where the primes denote derivatives with respect to u. This is a second order linear 
differential equation with regular singular points in the physical domain at the horizon 
u = 1 and the boundary u = 0. Solving equ. (11621) near the black hole horizon u — 1, 
we determine that the fluctuation of the metric is a linear combination of two solutions, 
h xy (co, u) ~ (1 — u) tiuj / 4:7tT . These solutions describe the gravitational wave propagating 
into (— ) and out of (+) the black hole, respectively. The infalling solution is the 
physically relevant retarded solution. Near the boundary u — > (or r — > oo) the 
gravitational field is also a linear combination of two solutions 

h xy (cu, u) = h° xy {uj) (1 + . . .) + B{uj) u 2 (1 + . . .) , (163) 

where . . . denotes terms that vanish as u — > 0. The two modes are called the non- 
normalizable mode and the normalizable mode. The non-normalizable mode is constant 
as r — > oo while the normalizable mode falls as 1/r 4 . Inserting the metric perturbation 
equ. (11631) into equ. ( 11521) the average stress tensor is 

(T xy (u)) =pK y {u) + (e + p) B(u) , (164) 

with the previously defined energy density and pressure, equ. ( 11531) . We observe that 
the coefficient of the non-normalizable mode, h° xyl can be interpreted as the external 
gravitational field applied to the gauge theory, while the coefficient of the normalizable 
mode, B(uj) is proportional to the induced stress tensor in the boundary theory. 

For an arbitrary value of B(uj) the general linear combination of solutions near the 
boundary would approach a linear combination of the infalling and outgoing solutions 
near the horizon. Thus the coefficient B(u) should be adjusted so that only the infalling 
solution (1 — u 2 )~W 47rT j s present near u — 1. In general the required B(to) has to be 
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Figure 5. Spectral function p xyxy (oj, k = 0) associated with the correlation function of 
the xy component of the energy momentum tensor. The spectral function is normalized 
to entropy density s. Left panel (Fig. (a)): Schematic picture of the spectral density 
in weak coupling QCD or SUSY Yang Mills theory [H51[FH)] . Right panel (Fig. (b)): 
Spectral density in strong coupling SUSY Yang-Mills theory calculated using the 
AdS/CFT correspondence, from |147j . 



determined numerically. For small u however, a straightforward calculation shows that 
to linear order in u the solution which is infalling at the horizon is 



h X y = K y {u) (1 - U)-^ T 



l-£^log(l + M ) + 0(^ 2 ) 



(165) 



Expanding this functional form near the boundary we find B(oS) = —iu}/(4irT). Then 
using e + p = sT and comparing the functional forms in equ. (11601) and equ. (j!64p we 
conclude that (T xy (u)) = ph° xy — iurjh° xy with 

- = ^-- (166) 

S 47T 

Remarkably, the strong coupling limit of the shear viscosity is small and independent 
of the coupling. The difference as compared to the weak coupling result becomes even 
clearer if one considers the spectral function. As described in Sect. I3.4l the Kubo formula 
relates the shear viscosity to the zero energy limit of the stress-energy spectral function. 
In weak coupling QCD the spectral function has a narrow peak near zero energy which 
reflects the fact that momentum transport is due to quasi-particles that are almost on- 
shell. The height of the transport peak is governed by the kinetic theory result for the 
shear viscosity. Kubo's formula implies that p(uj)/uj ~ T 3 /g 4 as uj — > 0. The width can 
be reconstructed from the /-sum rule 

r/ A ^ M = r <% +p) , (167) 

Jo u 5 
where g 4 T « A « g 2 T. Since the height of the transport peak is T 3 /g 4 , the width 
must be g 4 T. The high energy part of the spectral density can be computed from the 
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one-loop correlation function. The result is p(u) ~ uj a . A schematic picture of the 
spectral function is shown in Fig. [5](a). 

In the strong coupling limit the width of the transport peak becomes bigger, and 
the height becomes smaller. In M = 4 SUSY Yang Mills the infinite coupling limit 
can be determined as outlined above [1471 1148] . Specifically, the spectral function may 
be found by determining Gr(uj) from the numerical coefficient B(uj). The result is 
shown in Fig. [5](b). Clearly, the transport peak has completely disappeared, and there 
is no possibility of a quasi-particle interpretation of momentum transport. Whether 
the spectral function of the quark gluon plasma near T c looks more like Fig. [5^a) or 
(b) will have to be settled by numerical calculations on the lattice, see Sect. 13.41 It 
is interesting to note that numerical calculations of the shear viscosity, which require 
the determination of the zero energy limit of p(u)/uj, are easier in the case of strong 
coupling than they are for weak coupling. 



4.3. The KSS bound 

The calculation of the shear viscosity has been extended to other strongly coupled field 
theories with gravitational duals. It was discovered that within a large class of theories 
the strong coupling limit of r\ depends on the theory, but the ratio rj/s does not. This 
observation can be understood using Kubo's formula and the optical theorem. The 
optical theorem implies that the imaginary part of a correlation function can be related 
to the total cross section. As a consequence, the shear viscosity can be expressed in 
terms of the total graviton absorption cross section [19] , 

V (168) 

The low energy limit of cr a b s is equal to the area A of the event horizon, and the entropy 
density is given by the Hawking-Bekenstein formula, s = A/ (AG) where G = k|/(87t). 
The ratio rj/s is independent of A and k 5 . More formal arguments for the universality 
of r)/s in the strong coupling limit of field theories with holographic duals were given 
in |149[ [28l 1150] . Corrections to the infinite coupling limit of M = 4 SUSY Yang Mills 
theory were studied in [211 H51[ I152[ 1153] . The result is 

1503) 

A3/2 

The first correction is positive, as one would expect from the fact that rj/s — > oo as 
A — > 0. Based on these observations, Kovtun, Son, and Starinets (KSS) conjectured 
that 

- > ^ (170) 
s Air 

is a universal bound that applies to all fluids [10J. There is no proof of this conjecture, 
and a number of authors have attempted to construct counter examples. One possibility 
is a weakly interacting non-relativistic fluid with an exponentially large number of 
species or internal degrees of freedom, and therefore a very large entropy [TOl 11541 1155] . 




+ ...}. (169) 
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These systems are unusual because the time scale for thermal equilibration vastly 
exceeds the time scale for momentum equilibration, and because the fluid is not stable 
on very long time scales [156j . More recently, it was realized that theories with 
holographic duals described by higher derivative gravity may violate the KSS bound 
[1571 11581 11591 1160] . An explicit example was constructed by Kats and Petrov [159j . 
They showed that in Af = 2 SUSY Sp(N c ) gauge theory with a certain combination of 
matter fields 

^fl-^V (171) 



s 4tt V 2N C/ 

up to corrections of 0(A™ 3 / 2 ). For A 3 / 2 > iV c > A > 1 we find a violation of the 
KSS bound in a controlled calculation. However, there are bounds on the coefficients of 
higher derivative terms, and a modified bound on rj/s may yet exist [160J. 

4-4- Other transport properties 

There has been a large amount of work on applications of the AdS / CFT correspondence 
to transport properties other than the shear viscosity. Here we briefly summarize some 
results relevant to this review. M = 4 SUSY Yang Mills theory has a conserved R- 
charge (see Sect. 12. 5p . and we can study transport in the presence of a finite R-charge 
density. Son and Starinets find that the shear viscosity and entropy density depend 
on the density, but the ratio 77/s does not [161] . They also determine the thermal 
conductivity 

8tt 2 T 



/< 2 



V, (172) 



as well as the R-charge diffusion constant. The heavy quark diffusion constant was 
calculated in [1621 11631 1164j . The result is 

which depends on the value of the coupling A, and goes to zero in the strong coupling 
limit. The functional dependence on A is unusual from the point of view of perturbation 
theory, but typical of other AdS/CFT results. We also note that in the strong coupling 
limit the ratio of the heavy quark diffusion coefficient to the kinematic viscosity r]/(sT) 
goes to zero, whereas this ratio is independent of the coupling in the perturbative limit 
The bulk viscosity of M = 4 SUSY Yang Mills theory vanishes, but non-conformal 
deformations of the original AdS/CFT correspondence have been studied. Buchel 
proposed that in holographic models there is a lower bound on the bulk viscosity, 
( > 2(| — c 2 )?7, where c s is the speed of sound [165] . Note that the weak coupling 
formula involves the square of (~ — c 2 ). Gubser et al. considered a number of model 
geometries tuned to reproduce the QCD equation of state, and find that (/s has a 
maximum near the critical temperature where C/ s — 0-05 [166J. Larger values of £/s 
near T c have been suggested based on lattice data for the QCD trace anomaly [1£7_|. 
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4-5. Hydrodynamics and holography 

Up to this point we have used the AdS/CFT correspondence to calculate the transport 
coefficients that appear in first order hydrodynamics. However, AdS/CFT can be used 
to compute the full correlation function, and not just the hydrodynamic limit. An 
example is the spectral function shown in Fig. [51 and similar calculations have been 
performed in other channels as well. In this section we wish to discuss how the stress 
tensor of the fluid relaxes to the Navier-Stokes form. This process can be described by 
the second order terms introduced in Sect. 13. 1.31 We will follow the method outlined 
in |7Z|. 

A static fluid at temperature T corresponds to a black hole with a Hawking 
temperature Tjj = r /(7rL 2 ). First we switch to Eddington-Finkelstein coordinates, 
defining a new time coordinate v — t + J r dr L 2 /fr 2 . Then the metric is regular at the 
event horizon of the black hole, 

,.2 



r 



ds = 2 dv dr , , 



f(r)dv 2 + r 2 dx 2 . (174) 



Introducing four dimensional coordinates x^ = (x°, x 1 , x 2 , x 3 ) = (f,x), a vector 
= (1, 0) characterizing the local rest frame, and a scale parameter b characterizing 
the temperature the metric becomes 



^2 

ds 2 = -lu^dx^dr + — 

Lj 



-f{br)u^u v dx»dx v + r 2 P^ v dx l *dx v , (175) 



where P^ v = u^u v + rj^. The basic idea is to promote the variables u M and b to slowly 
varying functions of The metric is then 

ds 2 = — 2u fJ ,(x)dx' J 'dr 

y*2 



+ 15 



- f (b(x)r)u ll (x)u u (x)dx fl dx u + r P fJLV {x)dx^'dx v 

+ corrections due to gradients . (176) 

Variations in and b correspond to fluctuations in the local fluid velocity and 
temperature. Substituting this form into the Einstein equations, the corrections to 
the metric are determined order by order in the gradients of u^(x) and b(x). These 
metric corrections lead to deviations of the boundary stress tensor from an ideal fluid 
of precisely the form required by hydrodynamics. Up to second order we can write 

= Tq U + 5 W T» V + 5 {2) T^ + . . . , (177) 

and each term has physical significance. At zeroth order 

T o U = ^\^ T f W + ^u v ) , (178) 

which shows that e = 3P and that the pressure is 3/4 of the Stefan Boltzmann value. 
At first order 

N 2 

5 (l) T v» = __JL( 7rT )3 (7 ^ ) ( 179 ) 
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where o^ v is defined as in equ. (18"Uj) . This results shows that 77 = N 2 nT 3 /8. Combined 
with the zeroth order stress tensor we find rj/s — 1/47T, in agreement with previous 
results. Finally, at second order 



VTJI 



<D<r^ + -o^(d-u) 
3 



(180) 



where D = u ■ d, and the vorticity Vt^ as well as the transverse traceless tensor 

are defined in Sect. 13.1.31 The form of Tj% agrees with the general second order result 

for a conformal relativistic fluid derived in [75]. The second order coefficients are 

2-ln2 , 2ri , 2r?ln2 , , . 

r„ = — , Xl = JL, A 2 = ^, A 3 = 0. (181) 

We observe that the relaxation times are of order (7rT) _1 , the shortest time scale 
characterizing the plasma. 

4-6. Non-relativistic AdS/CFT correspondence 

Given the role that the AdS/CFT correspondence has played in improving our 
understanding of conformal relativistic fluids it is natural to ask whether the 
correspondence can be extended to non-relativistic scale invariant fluids like the dilute 
Fermi gas at unitarity. There has recently been significant progress in constructing 
holographic duals for non-relativistic field theories |168[ 11691 11701 1171j . 

The basic idea proposed in [1681 1169j can be explained by looking at the metric of 
d + 2 dimensional flat space 

ds 2 = Tj^dx^dx" = —2dx + dx~ + dx l dx l , (182) 

where we have introduced light cone coordinates ( ) with X ± = (x°±x d+1 )/V2 

and % = 1, . . . , d. Consider the massless Klein-Gordon equation in this space. In light 
cone coordinates 

(-2^ + gaW)^) = 0. (183) 

If the a;~-direction is compactified, then the corresponding momenta become discrete. 
We may write the lowest mode as 4>(x) ~ e -* ma; ^(x + ,Xj) and the equation for if) 
becomes the non-relativistic Schrodinger equation 

\2im 1 ^- + V 2 )il)(x + } x i ) = 0, (184) 



where x + plays the role of time. The symmetry group of this equation is known as the 
Schrodinger group Sch(d). The generators of the Schrodinger algebra include temporal 
and spatial translations, rotations, Galilean boosts, non-relativistic dilatations (which 
scale space and time by different factors, x — > sx and t —* s 2 t), a special conformal 
transformation (which scales t — ► t/(l + Xt) and x — > x/(l + At)), and the mass operator 
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The goal is to extend this construction to spaces that are asymptotically Anti- 
deSitter. The specific proposal in [1681 1169] is that the Schr(d) symmetry of a non- 
relativistic d + 1 dimensional conformal field theory can be mapped onto the isometries 
of the d + 3 dimensional metric 

ds 2 = r 2 (-2dx + dx~ - (3 2 r 2 {dx + ) 2 + (dx') 2 ) + (185) 

which reduces to the metric of AdSd+3 f° r /3 0. This metric can be realized in 
string theory by starting from AdSc, x X 5 , where X 5 is a generalized sphere called 
an Einstein-Sasaki manifold, and by applying a certain series of transformations that 
preserve solutions of the Einstein equations [1701 11711 1173j . The resulting field theory is 
a 2 + 1 dimensional field theory with infinitely many bosonic and fermionic fields, and 
an unusual equation of state P ~ T 4 //z 2 [173j . This is still quite far from the 3 + 1 
unitary Fermi gas, but the theory provides an explicit realization of a non-relativistic 
fluid which satisfies rj/s = l/(4-7r). 

The hydrodynamics of a holographic fluid with Schrddinger symmetry was studied 
in more detail in |174j . An interesting observation that was made in this paper is that 
the light cone reduction of a viscous relativistic stress tensor automatically leads to a 
VT term in the non-relativistic energy current. The thermal conductivity is completely 
fixed by the shear viscosity and the equation of state, 

k = 2 V 6 -±^ . (186) 

This result can be expressed in terms of the Prandtl number Pr = c p i]/ k, see equ. (I140p . 
Using the equation of state of a non-relativistic conformal fluid we find Pr = 1. The 
Prandtl number of many gases is indeed close to one, see Sect. 13.6.31 tmt at strong 
coupling there is no obvious reason for the relation Pr = 1 to hold. 



5. Experimental determination of transport properties 

In this section we will review experimental determinations of transport properties of 
liquid helium, cold atomic gases, and the quark gluon plasma. We will focus on shear 
viscosity, since it is the main focus of this review, and since it is the only transport 
property for which good data is available for all three systems. 

Liquid helium can be produced in bulk, and transport properties can be measured 
using methods that were developed for classical fluids. Cold atomic gases are produced 
in optical or magneto-optical traps. These traps typically contain 10 5 — 10 6 atoms. 
Hydrodynamic behavior is observed when the trapping potential is modified, or if the 
local density or energy density is modified using laser beams. The quark gluon plasma 
can only be created for brief periods in collisions of ultra-relativistic heavy ions. The 
system typically contains on the order of 10 3 — 10 4 quarks and gluons, and lasts for about 
10 fm/c (3 ■ 10~ 23 sec). Hydrodynamic behavior may take place during the expansion 
of the system and is reflected in the momentum spectra of particles in the final state. 
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Figure 6. Viscosity of 4 He at atmospheric pressure as a function of temperature. 
Data taken from Woods and Hollis-Hallett (green squares) [175 and Heikkila and 
Hollis-Hallett (red circles) [176j . This figure was adapted from [8]. The solid line 
shows the theory of Landau and Khalatnikov. The viscosity minimum corresponds to 
rj/n ~ 0.5 and rj/s ~ 1.9. 

5.1. Liquid helium 

There are a number of techniques for measuring the viscosity of fluids. Three popular 
instruments are: 

(i) Capillary viscometers are based on Poiseuille flow. Poiseuille's formula states 
that the flow through a pipe is inversely proportional to the shear viscosity, and 
proportional to the pressure drop as well as the fourth power of the diameter. 

(ii) Rotation viscometers measure the torque on a rotating cylinder or disk. The 
torque per unit length exerted by a pair of coaxial infinitely long cylinders is 
proportional to the shear viscosity and the difference between the angular velocities, 
and proportional to the ratio i?ii? 2 /(i?^ — R^), where i?i j2 are the two radii. 

(iii) Vibration viscometers determine the damping of an oscillating sphere or plate. 
These devices have many advantages but the data are more difficult to interpret, 
because the damping depends not only on the viscosity, but also on the density of 
the fluid. 

Initial measurements of the viscosity of superfluid liquid helium lead to an apparent 
contradiction between the results obtained using different methods. Capillary flow 
viscometers indicated vanishing viscosity below T c [177j . oscillatory viscometers showed 
a drop of the shear viscosity [178J, and experiments with rotation viscometers yielded 
a rise in viscosity below T c [175] . The contradictions can be resolved using superfluid 
hydrodynamics. The flow through a narrow capillary is entirely a superflow, and not 
sensitive to viscosity. Oscillation viscometers measure the product of viscosity and 
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normal density, which drops with temperature. Modern measurements confirm the rise 
of viscosity below T c which is predicted by the phonon-roton theory, see Fig. [61 The 
minimum viscosity of helium at normal pressure occurs just below the A point where 
i] ~ 1.2 • 10 -5 Poise. The minimum of rj/s occurs at higher temperature, close to the 
liquid gas phase transition. Recent measurements confirm the (weak) divergence of 
the shear viscosity at the critical endpoint of the liquid gas phase transition predicted 
by dynamical universality [179j . Experiments also find the expected (much stronger) 
divergence of the heat conductivity near the lambda point [H] . 

Once the shear viscosity and the heat conductivity are determined sound 
attenuation experiments can be used to measure the bulk viscosity (£2 in the superfluid 
phase) [180] . Below the A point (2 — 10~ 4 Poise. Damping of second sound determines 
a linear combination of (1 and (3 in the superfluid phase [181] , but the remaining linear 
combination is poorly constrained. 

5.2. Cold atomic gases 

Dilute Bose or Fermi gases are studied using optical traps that provide an approximately 
harmonic confinement potential 



The equilibrium density no can be determined from the equation of hydrostatic 
equilibrium, VPo — —n^V. Using the Gibbs-Duhem relation dP = ndfi + sdT we 
can see that this equation is solved by tiq(x) = n(fi(x)), where n(/z) is the equilibrium 
density as a function of the chemical potential, and fi(x) = fi — V(x). This result is 
known as the local density (or Thomas-Fermi) approximation, introduced by Thomas 
and Fermi in connection with the structure of heavy atoms. For dilute fermions at 
unitarity the equation of state at zero temperature is given by equ. (1271) and 



where \i is the chemical potential and 7 = 2/3. The chemical potential is related 
to the Fermi energy by the universal parameter £ introduced in equ. (12B1 . Transport 
properties of strongly interacting dilute Fermi gases can be extracted from a variety 
of experiments, free expansion from a deformed trap (elliptic flow) |182j . damping of 
collective oscillations [183[ 11841 [501 11851 [T86J, sound propagation |187j . and expansion 
out of rotating traps [188] . In the following we shall concentrate on damping of collective 
oscillations, as these experiments have been most carefully analyzed [1891 11901 11911 fl92j. 

We consider small oscillations around the equilibrium density, n = n + 5n. Since 
the damping is small, the motion is approximately described by ideal hydrodynamics. 
The compressibility at constant entropy is 






2 • 



(188) 




(189) 
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From the linearized continuity and Euler equation we get |193] 

m ^T = "7 (V • v) (W) -V(v-VV), (190) 

where we have dropped terms of the form VjVjV that involve higher derivatives of the 
velocity. This equation has simple scaling solutions of the form V{ = aiXi exp(iut) (no 
sum over i). Inserting this ansatz into equ. ( 11901) we get an equation that determines 
the eigenfrequencies u. The experiments are performed using a trapping potential with 
axial symmetry, oj\ = uj 2 = ujq, uj^ = Xu Q . In this case we find one solution with 
to 2 = 2u 2 and two solutions with [Ml EMI ESS] 

^ = ^ 2 { 7 + l + ^A 2 (191) 



In the limit of a very asymmetric trap (A — > 0) the eigenfrequencies are u 2 = 2uJq and 
uj 2 = (10/3)o;o- The mode u 2 = (10/3)cj 2 is a radial breathing mode with a = (a, a, 0) 
and the mode uj 2 = 2uJq corresponds to a radial quadrupole a = (a, —a, 0). 

The prediction of ideal hydrodynamics for the frequency of the radial breathing 
mode agrees very well with experimental results [183] . Damping of collective modes is 
due to viscous effects. The dissipated energy is given by 

E = - J d 3 x I (v iVj + y jV . - ls i:j V ■ v) 2 (192) 

+ CW(V-v) 2 + ^(VT)j, 

where r](x),((x) and k(x) are the local shear viscosity, bulk viscosity, and thermal 
conductivity. In the unitarity limit the system is scale invariant and the bulk viscosity 
in the normal phase vanishes. In the superffuid phase there are three bulk viscosities, 
Ci; C25 C3? see ec L u - flSD and (jT2j) . Scale invariance implies Ci = C2 = 0, see Sec. 13.6.11 and 
the contribution of £3 vanishes if v s = v n . For isentropic oscillations 5T ~ (5n/n)T. The 
solutions of equ. (I190p satisfy 8n(x) ~ n (x). This implies that there are no temperature 
gradients, and that thermal conductivity does not contribute to dissipation. 

We conclude that damping is dominated by shear viscosity. The energy dissipated 
by the radial scaling solutions is 



3 



(a 2 x + a 2 - a x a y ^ J d 3 xr](x), (193) 

where E is a time average. The damping rate is given by the ratio of the energy 
dissipated to the total energy of the collective mode. The kinetic energy is 

E kin = y J d 3 xn(x)v 2 = ^- (a 2 x + a 2 y ) (x 2 ). (194) 

In the case of a harmonic trapping potential the average (x 2 ) can be extracted using a 
virial theorem, E = 2N(V) [196] . The damping rate is 



IE 2 a 2 + a 2 — a x a y J d 3 xrj(x) 
2E ~ 3 a 2 + a 2 mN(x 2 ) 



(195) 
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Figure 7. Viscosity to entropy density ratio of a cold atomic gas in the unitarity 
limit, from [191] . This data points are based on the damping data published in |185j 
and the thermodynamic data in [1981 1199] . The light blue band is an estimate of the 
systematic uncertainty due to the breakdown of hydrodynamics near the surface of 
the cloud. The solid red and blue line show the low and high temperature limits of 
rj/s, see equ. (|123[) and (| 1 29[) . The dashed line shows the conjectured viscosity bound 
rj/s = l/(47r). 



where the factor 1/2 takes into account that the experiments measure an amplitude, not 
energy, damping rate. We note that the second factor on the RHS is 1/2 for the radial 
breathing mode and 3/2 for the radial quadrupole mode. This dependence provides an 
important check for the assumption that damping is dominated by shear viscosity. Also 
note that if the shear viscosity is proportional to the density or the entropy density then 
E/E scales as iV -1 / 3 . Near the surface the density is small and r}(x) will approach the 
Boltzmann limit, which is independent of density, see equ. (11291) . This is a problem, 
because the total volume of the system is infinite (at non-zero temperature the density 
has an infinite range tail). This difficulty is related to the breakdown of hydrodynamics 
near the surface of the cloud. An elegant solution to the problem is to include a finite 
relaxation time T v (r) = rj/(n(r)kBT) which diverges in the low density limit |197j . 

In order to compare with the proposed viscosity bound we will assume that the shear 
viscosity is proportional to the entropy density, rj(x) = as(x). Note that in general a 
is a function of T/Tp and varies across the trap (Tp depends on density). This means 
that we will extract an average value of a — rj/s. We can write 




e /2 (W) 1/S ^ ^ , (196) 



s 4 \u 2 

where T/oo± is the dimensionless damping rate, lo = ur^uj],^ is mean trap frequency, 
(S/N) is the entropy per particle, and E/Et=q is the equilibrium energy of the cloud in 
units of the zero temperature value. Fig. [7] shows rj/s extracted from the experimental 
results of the Duke group [185] . The entropy per particle was also taken from experiment 
[199J. Similar results are obtained if the entropy is extracted from quantum Monte Carlo 
data. We observe that rj/s in the vicinity of the transition temperature is about 1/2. We 
also note that the extracted shear viscosity roughly agrees with the high temperature, 
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Fermion quasi-particle, kinetic theory result. The low temperature, phonon dominated, 
result is not seen in the data, presumably because the phonon mean free path is bigger 
than the system size. 

There are many caveats that one should keep in mind regarding this analysis. 
First, we assume that shear viscosity is the only source of dissipation. There is 
some evidence for this assumption from comparisons of the damping rate of different 
collective modes [200] . On the other hand, the dependence of the damping rate on 
particle number predicted by viscous hydrodynamics has never been demonstrated. 
Second, hydrodynamics can only be applied in a relatively narrow temperature regime 
T < (2 — 3)T C . For higher temperatures the observed frequencies cross over from 
hydrodynamic behavior to a weakly collisional Boltzmann gas. This means that the 
kinetic theory prediction for the shear viscosity, equ. ( 11291) . is reliable but the frequency 
of the collective mode is too large for hydrodynamics to be applicable. Finally, there 
is an issue that is specific to the scaling flows (vi ~ a^Xj) considered here. Since the 
velocity field is linear in the coordinates, the second derivative of the velocity vanishes. 
This means that the viscous term in the Navier Stokes equation, pvi ~ Vj[i](ViVj + . . .)], 
is only sensitive to the density dependent part of the viscosity. But for a dilute gas the 
viscosity is expected to be density independent, see equ. (11291) . so the dilute limit can 
not be verified using experiments that involve scaling flows. 

There is clearly a need for additional experimental constraints. The first indication 
of almost ideal hydrodynamic behavior was the observation of elliptic flow by O'Hara 
et al. |182j . The experiment showed that if the trapping potential is removed the gas 
expands rapidly in the transverse direction while remaining nearly stationary in the 
axial. This is a consequence of the much larger pressure gradient in the short direction. 
The ideal hydrodynamics of this experiment was worked out in |201j but the effects 
of viscosity have not been carefully studied, in part because the data were taken at a 
single temperature. More recently Clancy et al. studied the expansion of a gas cloud 
with an initial velocity field corresponding to a scissors mode |188] . This is an interesting 
system, because the initial velocity field is irrotational (V x v = 0) but carries angular 
momentum. If the trapping potential is removed then the transverse size will grow 
initially, but if the gas remains irrotational then angular momentum conservation will 
force the transverse expansion to slow down (and the rotation to speed up) before the 
transverse and axial radii become equal [202] . This phenomenon was observed in the 
experiment, and an initial analysis leads to values of rj/s close to 1/(Att) [203, 204]. This 
result is very important, but some of the caveats mentioned above still apply. 

5.3. The quark gluon plasma at RHIC 

Cold quantum fluids can be studied in conditions that are very close to equilibrium. 
The quark gluon plasma, on the other hand, can only be created in relativistic heavy 
ion collisions. In these collisions the initial state is very far from equilibrium, and 
the system size is limited by the size of the heaviest stable nuclei. The applicability 
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of hydrodynamics is not clear a priori. In this section we will summarize some of 
the evidence that has been obtained from experiments at the Relativistic Heavy Ion 
Collider (RHIC). These experiments indicate that local equilibration takes place, that 
nearly ideal fluid dynamics is applicable, and that the shear viscosity to entropy density 
ratio near T c is within a factor of a few of the KSS bound. 

The collision energy in Au+Au collisions at RHIC is 100 GeV per nucleon, and the 
nuclei are Lorentz contracted by a factor of 7 ~ 100. The transverse radius of a Au 
nucleus is approximately 6 fm/c and on the order of 7000 particles are produced overall. 
The motion of the particles is relativistic, and the duration of a heavy ion event is r ~ 6 
fm. In order for hydrodynamics to be applicable this time has to be large compared to 
the equilibration time. 

The main observables are the spectra dN/ d 3 p of produced particles. For momenta 
less than 2 GeV the spectra roughly follow Boltzmann distributions with a characteristic 
temperature close to the QCD critical temperature. The first hint that the system is 
behaving collectively is the existence of radial flow. Heavy particle spectra have apparent 
temperatures that are larger than the temperatures extracted from light particles. This 
can be understood if there is a collective transverse expansion velocity v± which boosts 
the observed transverse momenta by an amount p± ~ mv±. 

More dramatic evidence for hydrodynamics is provided by the observation of elliptic 
flow in non-central heavy ion collisions. The centrality of the collision is characterized 
by the impact parameter b, the transverse separation of the two nuclei. The magnitude 
of b can be determined experimentally by selecting events with a given multiplicity of 
produced particles. The uncertainty in the impact parameter determination is small 
except in very peripheral bins [205] . The direction of the impact parameter can be 
determined on an event by event basis using the azimuthal dependence of the spectra. 
Imagine that the impact parameter direction is already known. This defines a coordinate 
system where z is along the beam axis and x is along the impact parameter direction, 
see Fig. [8j We write {p x ,P y ,Pz) = (p± cos(0),p_i_ sin(0),p 2 ), and the particle distribution 
can be expanded in Fourier components of 0, 



For a typical mid-central collision with b ~ 6 fm the v 2 harmonic, called the elliptic flow 
coefficient, is approximately 6%. In an actual event the reaction plane can be determined 
(in principle) by plotting the distribution in relative to an arbitrarily chosen axis, and 
then requiring that the distribution has a maximum at <ft = 0. This intuitive method 
to determine the reaction plane forms the basis of the event plane method. The result 
can be corrected for t> 2 fluctuations and additional correlations among the produced 
particles. Current analyses are not based on the event plane method but use two, 
four, and higher particle cumulants - see |206] and references therein for a complete 
review. The scaling of these cumulants with multiplicity demonstrates that one can 
reliably extract collective flow down to small system sizes. These measurements provide 



dN 



Vo(p ± )(l + 2v 2 ( P± ) cos(20) + 2v 4 (p ± ) cos(40) + ...). (197) 
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Figure 8. Geometry of a high energy heavy ion collision. The left panel shows the 
collision of two Lorentz contracted gold nuclei. The beam direction is the z-axis. The 
right panel shows the same collision in the transverse plane. The impact parameter is 
along the x-axis, and the remaining transverse direction is the y-axis. 

a unique opportunity to study the approach to hydrodynamic behavior in a controlled 
fashion [207] . 

Elliptic flow represents the collective response of the system to pressure gradients 
in the initial state. At finite impact parameter the initial state has the shape of an 
ellipse, with the short axis along the x direction, and the long axis along the y direction. 
This implies that pressure gradients along the x axis are larger than along the y axis. 
Hydrodynamic evolution converts the initial pressure gradients to velocity gradients in 
the final state. Elliptic flow is a direct measure of collectivity. In particular, if the 
nucleus nucleus event were a simple superposition of proton proton collisions then the 
particle distribution would be azimuthally symmetric. 



5.3.1. The Bjorken model The application of hydrodynamics to relativistic heavy ion 
collisions goes back to the work of Landau |208] and Bjorken [209J. Bjorken discovered a 
simple scaling solution that provides a natural starting point for more elaborate solutions 
in the ultra-relativistic domain. Consider two highly relativistic nuclei moving with 
equal but opposite momenta in the z direction. In the relativistic regime the natural 
variable to describe the motion in the z direction is the rapidity 

V=k**(¥*)- (198) 



E-p, 

At RHIC the energy of the colliding nuclei is 100+100 GeV per nucleon, and the 
separation in rapidity is Ay = 10.6. Bjorken suggested that the two highly Lorentz 
contracted nuclei pass through each other and create a longitudinally expanding fireball 
in which particles are produced. In the original model the number of produced particles 
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is independent of rapidity, and the subsequent evolution is invariant under boosts along 
the z axis. The evolution in proper time is the same for all comoving observers. The 
flow velocity is 

Up = 7 (1, 0, 0, v z ) = {t/r, 0, 0, z/r), (199) 



where 7 is the boost factor and r = ^t 2 — z 2 is the proper time. The velocity 

field (11991) solves the relativistic Euler equation (178j) . In particular, there is no 

longitudinal acceleration. The remaining hydrodynamic variables are determined by 
entropy conservation. Equ. ( 1771) gives 

jp Mr)] = (200) 

and s(t) = s t /t. For an ideal relativistic gas s ~ T 3 and T ~ 1/r 1 / 3 . Typical 
parameters at RHIC are r D ~ (0.6 — 1.6) fm and T Q ~ (300 — 425) MeV. The combination 
t T 3 is constrained by the final multiplicity, but individually r G and T Q are not well 
constrained. We note that the corresponding initial temperature is significantly larger 
than the critical temperature for the QCD phase transition. 

The temperature drops as a function of r and eventually the system becomes to 
dilute for the hydrodynamic evolution to make sense. At this point, the hydrodynamic 
description is matched to kinetic theory, 

T ky d ro = T «n = J dT p ^ /(X) p t) (201) 

For an ideal fluid the distribution function is parameterized by the local temperature 
and flow velocity 

~ t exp(p ■ u/T) ± 1 

where i labels different particle species, and di are the corresponding degeneracies. 
Finally, the observed particle spectra are given by 

= ^-p J d£^(x,p,t), (203) 

where is the normal vector on the "freezeout" hypersurface, the surface on which 
the matching between the hydrodynamics and kinetic descriptions is performed. 

In order to quantitatively describe the observed particle distributions several 
improvements of the simple Bjorken model are necessary. First, one has to include 
the transverse expansion of the system [210] . Transverse expansion becomes important 
at a proper time r ~ R/c s , where R is the (rms) size of the nucleus, and c s is the speed 
of sound. At very late times the expansion becomes three dimensional, 

sir) ~ ^ > ( 204 ) 

and T ~ 1/r. Transverse expansion is caused by transverse pressure gradients. These 
gradients are sensitive to the initial energy density of the system. One simple model 
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for the initial energy density (or entropy density) in the transverse plane is the Glauber 
model. In the Glauber model the entropy density is 

s(x ± , b) oc T A (x± + b/2) 1 - exp (-o NN T A (x± - b/2)) 

+ T A (x ± - b/2) [l - exp (-a NN T A (x ± + b/2)) j , (205) 

where b is the impact parameter 

Ta(x_l) = J dzp A (x) (206) 

is the thickness function, and o~nn(Vs) is the nucleon-nucleon cross section. Here, Pa( x ) 
is the nuclear density. The idea behind the Glauber model is that the initial entropy 
density is proportional to the number of nucleons per unit area which actually collide. 
Other variants exist. For instance, one can distribute the energy density according to 
the number of binary nucleon-nucleon collisions - see Ref. [211] for a comparison. A 
more sophisticated theory of the initial energy density is provided by the Color Glass 
Condensate (CGC) [2121 1213] . This model leads to somewhat steeper initial transverse 
energy density distributions. 

Gradients in the transverse pressure lead to transverse acceleration and generate 
collective transverse flow. The collective expansion leads to a blue-shift of the transverse 
momentum spectra of produced particles. For an azimuthally symmetric source with 
temperature T and radial flow velocity u r (r) we get |214] 

dN\ 2di r fm^u 1 



Po 



d 3 pj, (2tt) 



fm ,Jrdr Kl (^f-) h (^f)^ (207) 

where m± = \Jp\ + mf is the transverse "mass", and Kx(z), I\{z) are generalized Bessel 
functions. Using the asymptotic form of the Bessel functions one can show that the 
spectrum has the form 



exp -— , T eff = TJ- . (208) 



dm{ 1 \ T eff J ' * JJ V l-«r 
This effect of transverse flow on the spectra is seen in the data. At RHIC, transverse 
velocities at freezeout reach 0.6 c. At finite impact parameter the initial energy density 
in the transverse plane is not azimuthally symmetric. The pressure gradient along the 
direction of the impact parameter is larger than the gradient in the orthogonal direction. 
The resulting anisotropy of the transverse flow can be characterized by the elliptic flow 
parameter t> 2 defined in equ. (j!97p . The observed elliptic flow is remarkable because t>2 
is a very direct measure of transverse pressure. The observed radial flow is proportional 
to the radial pressure and the expansion time, but elliptic flow has to be generated early, 
when the system is still deformed. 



5.3.2. Estimates of viscous corrections We are now in a position to discuss the role of 
dissipative effects. We begin with the effect of shear and bulk viscosity on the Bjorken 
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solution. The scaling flow given in equ. (11991) is a solution of the relativistic Navier- 
Stokes equation. Viscosity only modifies the entropy equation. We get 



Ids 1/ fr? + C 



(209) 



s dr t \ sTt 

We observe that dissipation is governed by the sound attenuation length T s , see equ. (151)1 . 
The applicability of the Navier-Stokes equation requires that the viscous correction is 
small [9] 

- + 7-«!(?Y)- (210) 
s 4 s 4 

For the Bjorken solution Tr ~ r 2 / 3 grows with time, and this condition is most restrictive 
during the initial phase. Using r D = 1 fm and T = 300 MeV gives 77/s < 0.6. For a 
three dimensional expansion Tr is independent of time. At very late time the fluid is 
composed of hadrons, or pre-formed hadronic resonances. In that case rj ~ T/a, where 
a is a hadronic cross section. Then, for a three dimensional expansion, the viscous 
correction r)/(sTr) grows with proper time as r 2 . This shows that the system has to 
freeze out at late time. 

It is instructive to study the viscous contribution to the stress tensor in more detail. 
At central rapidity we have (for ( = 0) 

T ZZ P ~ , T XX Tyy P + — . (211) 

O T O T 

This means that shear viscosity decreases the longitudinal pressure and increases the 
transverse one. In the Bjorken scenario there is no acceleration, but if pressure 
gradients are taken into account shear viscosity will tend to increase transverse flow. A 
similar effect will occur at finite impact parameter, see Fig. [8j Shear viscosity reduces 
the pressure along the x-direction, and increases the pressure in ^-direction. As a 
consequence there is less acceleration in the ^-direction, and elliptic flow is suppressed. 
This is the basic observation that motivates attempts to extract shear viscosity from 
the observed elliptic flow. 

Viscosity modifies the stress tensor, and via the matching condition (12011) this 
modification changes the distribution functions of produced particles. In Ref. |215] a 
simple quadratic ansatz for the leading correction 5f to the distribution function was 
proposed 

Sf = |^/o(l + /o)j>oP/jV<V> , (212) 

where V^u^ is a symmetric traceless tensor, see equ. (1861) . This form summarizes the 
results of more involved kinetic calculations [17J. The modified distribution function 
leads to a modification of the single particle spectrum. For a simple Bjorken expansion 
and at large p± we find 

sm _ r. fP.V (213) 



dN 4r/ V T 
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where ry is the freezeout time. We observe that the dissipative correction to the 
spectrum is controlled by the same parameter Y s /r that appeared in equ. (12091) . We also 
note that the viscous term grows with p±. At RHIC transverse momentum spectra are 
in agreement with hydrodynamic predictions out to transverse momenta several times 
larger than the temperature. In equ. (12131) this is partially compensated by the fact 
that Tf/r is a large number, but typically the requirement S(dN)/(dN ) < 1 provides a 
more stringent bound on rj/s than equ. (12101) . We can also include transverse expansion 
and study the leading dissipative correction to t>2 [215j . The viscous correction tends to 
reduce v 2 and grows with p±. At p± = 1 GeV an estimate similar to equ. (12 13[) gives 
(Sv2)/v2 ~ 1 for r s /r/ ~ 0.2. Using r/ ~ 5 fm and Tf = 160 MeV this translates into 
r)/s < 0.6. 

5.3.3. Hydrodynamic Simulations There have been a number of recent numerical 
studies devoted to extracting the shear viscosity of the quark gluon plasma |2161 12171 
12181 12191 1220] . Here, we will follow the work of Dusling and Teaney [216], and refer 
the reader to the recent review by Heinz [30] for a more detailed comparison between 
different strategies for implementing relativistic viscous hydrodynamics for heavy ion 
collisions at RHIC. In order to respect causality in viscous hydrodynamics we have to use 
a second order formalism. This means that the strains 8T^ U are promoted to dynamical 
fields which relax on a collisional time scale to the Navier-Stokes form rather than 
being specified instantaneously by the constitutive relations. Such relaxation processes 
are second order in the hydrodynamic expansion, see Sec. 13.1.31 Dusling and Teaney 
considered a 2 + 1 dimensional boost invariant flow and used a second order fluid model 
studied by Ottinger |221j . This model is formulated along the lines of equ. (1551) . and 
the dynamical strains are parametrized by an additional field (f v , 

ST 1 " = -P{e)acT. (214) 

The relaxation equation for (f ,v is written in terms of a tensor c^ u . This tensor satisfies 
the constraint 

c ilv u v = Up , (215) 
and is decomposed as 

Cfj,„ = - u^Uy + c^ + Cfu, , (216) 
- (c x x -l){ %lu + u,u u ) . (217) 
The field c^ v obeys the relaxation equation 

u x (d x c^ - d^c Xu - dyCfo) = -—Zfa, - —Cfa,. (218) 

The constraint equ. (12151) is preserved under time evolution. In the limit that the 
relaxation times tq and T2 are small the evolution equation leads to the following solution 
for c^ u in the local rest frame 

c « = r 2 (W + &u< - ^5 tj d k u k ^j + ~ r 6 ij d k u k . (219) 
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Figure 9. Plot of the energy density per unit rapidity er (left) and of the transverse 
velocity (right) at times of r = 1,3,6,9 fm/c, for rj/s = 0.2 (solid red line) and for 
ideal hydrodynamics (dotted blue line), from |216j . 



Comparing with the Navier-Stokes equation we see that 

r] = T 2 Pa, ( = ^r Pa. (220) 

Dusling and Teaney used these relations (with a = 0.7) to fix (r , r 2 ) in terms of r\ and 
C, and studied the sensitivity to the parameter a. They considered a simple conformal 
equation of state P = e/3, and set ( = 0. 

An advantage of Ottinger's approach is that equ. (12181) is relatively simple to solve. 
One just evolves the spatial components of Cjj and then uses the constraints (121 5[) to solve 
for the time components Coo and Coj. As hydrodynamics is universal, any fluid model 
can be recast in terms of the first and second order formalism described in Sec. 13.1.31 
For the Ottinger fluid model, expanding out the equations of motion to second order 
leads to the relation 



tT = - W ^ - T 2 



(221) 



+ ZW>A + r 2 7r ( ^> A - \ r 2 ^(d ■ u) 



The last term in this expression differs from the general result in equ. (I22ip . which 
indicates that at second order in gradients this model contains terms that break 
conformal symmetry. 

The hydrodynamic equations were solved for several fixed values of the shear 
viscosity to entropy density ratio rj/s. At the initial time r — 1 fm the entropy 
per participant is adjusted and closely corresponds to the results of full hydrodynamic 
simulations [2221 12231 [224J. The maximum initial temperature is T D = 420 MeV at an 
impact parameter 6 = 0. The initial components of the stress tensor are set to the 
Navier Stokes values. 

Numerical results are shown in Fig. [9j The effect of viscosity is twofold. The 
longitudinal pressure is initially reduced and the viscous case does less longitudinal PdV 
work as in the simple Bjorken expansion. This means that at early times the energy 
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Figure 10. Location of freezeout surfaces for for central Au-Au collisions. The surfaces 
are determined by the condition % (d ■ u) = 0.6, slightly larger than the value 0.5 
discussed in the text. Different surfaces correspond to different values of rj/s. The 
shading corresponds to the freezeout temperature. The thin solid black curve shows 
the contour set by ^ (d ■ u) — 0.225 for comparison. 

per rapidity decreases more slowly in the viscous case. The reduction of longitudinal 
pressure is accompanied by a larger transverse pressure. This causes the transverse 
velocity to grow more rapidly. The larger transverse velocity causes the energy density 
to deplete faster at late times in the viscous case. The net result is that a finite viscosity, 
even as large as rj/s = 0.2, does not integrate to give major deviations from the ideal 
equations of motion. 

Freezeout occurs when the viscous terms become large compared to the ideal terms. 
Roughly, the system begins to break apart when 

^-d-u~l. (222) 
P v ' 

This combination of parameters can be motivated from kinetic theory. The pressure 

is of order P ~ e (v^ h ) where (v% h ) is the typical quasi particle velocity and e is the 

energy density. The viscosity is of order rj ~ e (v^ h ) tr where r# is the relaxation time. 

Thus the freezeout condition is (d ■ u) ~ r R (d ■ u) ~ 1. Ideally there should be an 

overlap region where both viscous hydrodynamics and kinetic theory are valid. In this 

region hydrodynamics can be systematically coupled to a kinetic description in order 

to correctly model the breakup. In simulations of heavy ion collisions the size of the 

overlap region is small, and the breakdown of hydrodynamics is typically modeled via a 

freezeout surface. In practice the freezeout surface was chosen to satisfy (rj/P) d-u = 0.5, 

where the precise number on the right hand side is simply an educated guess based on 

examining the output of second order hydrodynamic simulations. Typical freezeout 

surfaces which satisfy this criterion are shown in Fig. [101 We note that hydrodynamics 
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Figure 11. Left: V2(pt) for massless Bose particles for simulations using rj/s = 
0, 0.05, 0.2 at an impact parameter of b = 6.5 fm . Right: Four-particle cumulant data 
as measured in Au-Au collisions at s/s = 200 GeV for a centrality selection of 16% to 
24% In- 



breaks down both at late and also at early times. The latter is most clearly see from the 
rj/s = 0.4 curve. If rj/s is very small then freezeout will occur at very late proper time. 
In the following, we shall therefore use a simpler criterion x= h(d ■ u) = const which 
is independent of rj/s. Taking x — 3 roughly corresponds to the rj/s = 0.2 freezeout 
surface in Fig. [TUJ 

Finally, we can compute the spectra of produced particles. We follow the procedure 
outlined above, see equ. (I212p . and write the distribution function as / = f + 5f with 

/ (1t/o)«,. 



Sf 



(223) 



2(e + P)T 2 

The spectrum is determined by integrating the distribution function over the freezeout 
surface as in equ. (12031) . and the elliptic flow parameter v 2 is computed from the definition 
in equ. (1197j) . A comparison with the data obtained by the STAR collaboration is shown 
in Fig. (ITT1) . There are several curves here and we will go through them one by one: 

• For the two different values of rj/s, 0.05 and 0.2, there are three curves each. Our 
best estimates for the elliptic flow as a function of pr are labeled f + 5f n and are 
shown by the blue (rj/s = 0.05) and orange (rj/s = 0.2) lines. 

• To disentangle what part of the viscous modification is due to the distribution 
function 5f and what part of it is due to changes in the flow, we also compute 
V2(pt) with /o only. We see that the effect on v 2 from viscous modifications of the 
flow is relatively minor. This may be the largest obstacle to reliably extracting the 
shear viscosity from the heavy ion data. However, it is important to realize that 
the modification to the distribution function reflects the viscous correction to the 
stress tensor itself. In hydrodynamic simulations the pt integrated v 2 tracks the 
asymmetry of the stress tensor |226j . In the context of viscous hydrodynamics this 
is nicely illustrated in Fig. 8 of Ref. [219J. Nevertheless, in contrast to the atomic 
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physics experiments discussed in Sect. I5.2[ the observed viscous corrections to the 
elliptic flow do not reflect a resummation of secular terms in the gradient expansion. 

• Finally instead of showing the spectrum computed with <5T Ml/ , we show v 2 (pt) 
computed with velocity gradients directly. For this purpose 5T^ U in equ. (12231) is 
replaced by —r/a^ where o^ v is computed from the flow velocities, v 2(pt) computed 
in this way is denoted by f + 5 fa in the figure and corresponds to the magenta 
and black curves. To first order in gradients this is an identity and the difference 
between f n and fa is a measure of the magnitude of second order terms. For the 
smallest viscosity rj/s = 0.05 the differences are quite small, but the effect is more 
noticeable for rj/s = 0.2. 

We conclude that for small values of rj/s < 0.2 the gradient expansion is working. 
There are, however, a number of issues that have to be considered in order to extract 
reliable values for r)/s: 

• The constraints on rj/s are sensitive to the initial values for the transverse energy 
density. In particular, using color glass initial conditions produces higher transverse 
pressure gradients, and allows for values of rj/s about twice larger than Glauber 
model initial conditions |219j . 

• Near the edge of the nucleus the gradient expansion breaks down completely. It is 
important to quantify the extent to which the effects of the edge propagate into 
the interior and invalidate the hydrodynamic description. One way to do this is by 
comparing the results of hydrodynamic simulations to kinetic theory. For strongly 
coupled plasmas kinetic theory is not an appropriate description of the microscopic 
interactions. However, hydrodynamics is independent of the microscopic details. 
Thus extrapolating kinetic theory into the strongly coupled domain is good way 
to construct a model which gracefully transitions from a hydrodynamic description 
in the interior to a kinetic description near the edge. There are several important 
developments in this direction [2201 12271 1228] . 

• Viscous effects in the hadronic phase are very important |222[ 1213] . These effects 
can be taken into account by coupling the hydrodynamic evolution to a hadronic 
cascade. 

• Effects due to bulk viscosity may reduce both radial and elliptic flow [2291 230J. 
Bulk viscosity is likely to be much smaller than shear viscosity in the plasma phase, 
but it is expected to grow near the phase transition. The magnitude of this growth 
is not clear. 

Some of these effects have yet to be carefully studied. However, even if one 
conservatively assumes that all uncertainties tend to increase the bound on rj/s, one still 
has to conclude that for a shear viscosity of rj/s > 0.4 it will be impossible to reproduce 
the observed flow. The question now is whether it will be possible to describe the large 
set of available data on energy, impact parameter, rapidity, transverse momentum, and 
species dependence of flow using viscous hydro, and whether it is possible to extract a 
reliable value, with controlled error bars, for rj/s of the quark gluon plasma. 
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6. Summary and outlook 

6.1. Summary 

In this review we summarized theoretical and experimental information on the behavior 
of nearly perfect fluids. We characterized the "perfectness" in terms of the shear viscosity 
of the fluid. Shear viscosity is special because 

• shear viscosity is the "minimal" transport property of a fluid. Thermal transport, 
diffusion, or conductivities require the presence of conserved charges, and bulk 
viscosity vanishes in the case of scale invariant fluids. The "perfectness" of any 
fluid can be characterized by the dimensionless ratio rj/s. If factors of h and ks are 
reinstated the dimensionless measure of fluidity is [rj/h]/[s/kB\. 

• a small shear viscosity is uniquely associated with strong interactions. Other 
transport coefficients may vanish even if the interaction strength remains weak. 
For example, bulk viscosity vanishes if the theory is conformal, and the diffusion 
constant goes to zero at the a localization phase transition, but the shear viscosity 
of a weakly coupled fluid is always large. 

Note that the reverse of the last statement is not true: The viscosity can be large, 
even if the interaction is strong. One possible scenario is that the interaction is so strong 
that it leads to the breakdown of a continuous symmetry, and the emergence of a new set 
of weakly coupled quasi-particles. Note that the M = 4 super-conformal fluid is special 
in this regard: There is no symmetry breaking in the strong coupling limit. Another 
example is the liquid-gas phase transition. There are no weakly coupled quasi-particles, 
but the viscosity diverges because of critical fluctuations. 

We summarized the main theoretical approaches to transport coefficients: kinetic 
theory, holography, and non-perturbative approaches based on the Kubo formula. 

• Kinetic theory applies whenever the fluid can be described in terms of quasi- 
particles. For many fluids this is the case in both the high and the low temperature 
limit. Typically, at high temperature the quasi-particles are the "fundamental" 
degrees of freedom (quarks, gluons, atoms) whereas the low temperature degrees of 
freedom are composite (phonons, rotons, pions). In the regime in which kinetic 
theory applies the ratio rj/s is always parametrically large. This leads to a 
characteristic "concave" temperature dependence of rj/s. Kinetic theory is useful 
in constraining the location of the viscosity minimum (usually, near the crossover 
between the high and low temperature regimes). Also, despite the weak coupling 
restriction, kinetic theory is quantitatively accurate for many quantum fluids in the 
whole temperature range covered by experiment. 

• Holography is a new method for studying the transport behavior of quantum 
fluids. It is most useful in the strong coupling limit of certain model field theories 
(like M = 4 super-conformal field theory), but the range of field theories that 
have known holographic duals has grown significantly over the years. More 
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importantly, holographic dualities have led to important new insights into the 
transport properties of strongly coupled fluids. This includes the proposed universal 
bound on rj/s, consistent higher order hydrodynamic theories, computations of 
the spectral function associated with the shear and other transport modes, etc. 
Holography has also been used to study specific solutions to the hydrodynamic 
equations, like the wake of a moving heavy quark, or the approach to equilibrium 
after the collision of two highly Lorentz contracted sources. 

• The Kubo formula connects non-perturbative calculations of equilibrium correlation 
functions with non-equilibrium transport coefficients. Possible non-perturbative 
approaches include Euclidean lattice, large-N, exact renormalization group, and 
Dyson-Schwinger calculations. Significant progress has been made in computing the 
shear viscosity of the quark gluon plasma on the lattice. The results are close to the 
proposed bound. Future calculations will answer the question whether transport 
phenomena in the QCD plasma can be understood in terms of quasi-particles. 

Finally, we summarized the experimental situation for the three most strongly 
coupled fluids that can be prepared in the laboratory. 

• Liquid helium has been studied for many years and its shear viscosity is well 
determined. The minimum value of rj/s is about 0.8 and is attained near the 
endpoint of the liquid gas phase transition. The ratio rj/n has a minimum closer 
to the lambda transition. Even though rj/s < 1 the transport properties of liquid 
helium can be quantitatively understood using kinetic theory. 

• Strongly interacting cold atomic Fermi gases were first created in the laboratory in 
1999. These systems are interesting because the interaction between the atoms can 
be controlled, and a large set of hydrodynamic flows (collective oscillations, elliptic 
flow, rotating systems) can be studied. Current experiments involve 10 5 — 10 6 atoms, 
and the range of temperatures and interaction strengths over which hydrodynamic 
behavior can be observed is not large. There are also some difficulties in extracting 
the viscosity that are related to the nature of the flow profiles that have been 
studied. A conservative estimate is rj/s < 0.5. 

• The quark gluon plasma has been studied in heavy ion collisions at a number of 
facilities, AGS (Brookhaven), SPS (CERN), RHIC (Brookhaven). Almost ideal 
hydrodynamic behavior was observed for the first time in 200 GeV per nucleon (in 
the center of mass) Au on Au collisions at RHIC. These experiments are difficult 
to analyze - the initial state is very far from equilibrium and not completely 
understood, final state interactions are important, and the size and lifetime of 
the system are not very large. Important progress has nevertheless been made in 
extracting constraints on the transport properties of the quark gluon plasma. A 
conservative bound is rj/s < 0.4, but the the value of rj/s that provides the best fit 
to the data is smaller, rj/s ~ 0.1. 
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6.2. Outlook 

Much work remains to be done in order to advance theoretical methods for predicting the 
transport properties of strongly coupled quantum fluids, to improve the determination 
of transport coefficients of these fluids, and to discover new nearly perfect fluids. 

However, in addition to that, we also want to understand what nearly perfect are 
fluids are like, in particular whether they can be understood using quasi-particles and 
the tools of kinetic theory. There are several avenues for addressing this question: 

• Quantum Monte Carlo calculations can be used to determine the spectral function 
of the energy momentum tensor. If energy and momentum are carried by quasi- 
particles then the spectral function has a peak at low energy. If, on the other 
hand, the fluid is AdS/CFT-like then there are no quasi-particles and no peak 
in the spectral functions. Current calculations in lattice QCD seem to prefer the 
AdS/CFT picture [97l 1103] . but the issue is far from settled. 

• Detailed simulations of kinetic equations can be used to study the crossover from 
kinetic behavior (the Knudsen limit) to hydrodynamic behavior, see [ 220 ] 12271 1231[ 
1228] for work in the context of QCD and |115U197j for studies of the dilute Fermi gas 
at unitarity. The main question is whether it is possible to extend a self-consistent 
kinetic theory into the domain rj/s < 1, and whether one can describe not only flow 
properties, but also all other transport properties like diffusion and energy loss. Self- 
consistency requires that the life time of the quasi-particles is long compared to the 
characteristic thermal time, r ~ 1/T. A simultaneous description of momentum 
diffusion, charge diffusion, and energy loss is important because in kinetic theory 
there is a close connection between these observables. In AdS/CFT-like fluids 
the relation between different transport processes is non-trivial. For example, at 
infinite coupling shear viscosity is finite while the heavy quark diffusion coefficient 
is zero. The current experimental situation in QCD is not entirely clear, with both 
AdS/CFT based [232] and kinetic approaches [233] demonstrating some success. 

• Phenomenological studies address the quasi-particle structure of nearly perfect 
fluids by studying fluctuations and correlations of conserved charges other than 
energy and momentum. In liquid helium quasi-particles were studied using neutron 
scattering, and in the dilute Fermi gas one can study the quasi-particle structure 
using radio- frequency spectra, see [2341 1235] . In the quark gluon plasma there 
is some evidence for quasi-particles from charge fluctuations [236] and from the 
success of the recombination model in reproducing the flow of identified particles 
at intermediate momenta [237J. 
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